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ABSTRACT 

We present Spitzer Infrared Spectrograph spectra of 28 Class I protostars in the Taurus star-forming 
region. The 5 to 36 /im spectra reveal excess emission from the inner regions of the envelope and 
accretion disk surrounding these predecessors of low-mass stars, as well as absorption features due to 
silicates and ices. Together with shorter- and longer-wavelength data from the literature, we construct 
spectral energy distributions and fit envelope models to 22 protostars of our sample, most of which are 
well-constrained due to the availability of the IRS spectra. We infer that the envelopes of the Class I 
objects in our sample cover a wide range in parameter space, particularly in density and centrifugal 
radius, implying different initial conditions for the collapse of protostellar cores. 
Subject headings: circumstellar matter — stars: formation — stars: pre-main sequence — infrared: 
stars 



1. INTRODUCTION 

In a now widely accepted evolutionary sequence based 
on the shape of the spectral energy distribution (SED) 
in the infrared, there are four stages in the pre-main- 
sequence life of a low-mass star: Class objects are 
deeply embedded protostars with large submillimeter to 
bolometric luminosity ratios, which are very faint at 
near-infrared wavelengths; Class I objects are surrounded 
by infalling envelopes, but are less embedded; Class II ob- 
jects have cleared their envelopes and are surrounded by 
accretion disks; and Class III objects h ave dispersed al- 
most all of their circumstellar mater ial (|Lada fe W ilking 
[1981: lLadal[T987l : lAndre et al.lll993h . 

Therefore, Class I protostars are thought to be in one of 
the earliest evolutionary stages. After the gravitational 
collapse of a molecular cloud core starts, the infalling ma- 
terial forms an envelope around the central object. Since 
cloud cores are usually rotating slowly, the material falls 
preferentially onto a disk, which serves as a mass reser- 
voir, most of which is eventually accreted by the proto- 
star. Class I objects are thus surrounded by accretion 
disks, in addition to envelopes extending to several 1000 
AU. The radius of the disk is roughly equal to the cen- 
trifugal radius (i? c ), which defines the centrifugal barrier 
encountered by matter falling in on the equatorial plane, 
if the specific angular momentum is conserved; inside this 
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radius the density distribution flattens, and the matter 
follows non-radial trajectories, falling onto a disk. Ma- 
terial tends to accumulate close to R c , which is close to 
the protostar at the begi nning of the collapse, but in- 
creases rapidly over time (jTerebev. Shu, fc CassenlllQSl 
I Adams fc Shul[l986h . 

Besides accreting matter, young stars generate pow- 
erful outflows which are launched along magnetic field 
lines; mass accretion onto the star and mass loss 
in the form of outflows seem to be correlated (e.g., 
iHartiganTE dwards . fc Ghandourl I1995T ) . In addition to 
collimated outflows, protostars might also emit a wide- 
angle wind, which is expected to op en up cavities in 
the envelope (|Arce fc Goodmaril 2002). The presence of 
both collimated outflows and wide-angle stellar winds ex- 
plains the fact that some objects with large cavities have 
well-collimated outflows. A partially evacuated cavity 
also f orms as a natural by-product of the collapse of a 
sheet (|Hartmann et al.lll994t lHartmann. Calvet. fc Boss! 
fl996h : since molecular cloud cores are rarely spherically 
symmetric, this scenario might be more realistic. In some 
cases, both stellar winds and the flattening of the enve- 
lope are responsible for cavities oriented along the rota- 
tional axis of a protostar. 

The spectral index n (XF\ cx A") from 2 to 25 fxm is 
used for the SED cl assification of protostars and pre- 
main - sequence stars (iLadal Il987t lAdams. Lada, fc Shu! 
Il987t I Andre fc Montmerlel Il994fk Class I objects have 
n > 0, Class II objects have — 2 < n < 0, and Class 
III objects have — 3 < n < — 2 (n = — 3 is the spec- 
tral index of a photosphere in the Rayleigh- Jeans limit). 
Therefore, Class I objects are characterized by a strong 
infrared excess; the dust in the envelope absorbs radia- 
tion from the central protostar and disk and reemits it 
at infrared to mm wavelengths. The peak of the SED 
of this type of objects usually lies in the mid- to far- 
infrared, depending on the density of the envelope (and 
thus m ass and infall rate), as well as the s ource lumi- 
nosity (|Kenvon. Calvet. fc Hartmannlll993"ah . Flat spec- 
trum sources are likely objects surrounded by envelopes 
and disks and therefore protostars, but with wide cavi- 
ties in their envelopes and/or seen close to pole-on (e.g., 
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DG Tau, T Taii; iFurlan et all [20061 : iCalvet et al.lfl99l . 
Some of the flat-spectrum sources could be transition ob- 
jects (Class I/II), consisting of a T Tauri star and its disk 
surrounded by remnant envelope material (e.g., IRAS 
04154+2823). Their SEDs typically show a double peak, 
one centered in the near-IR and one in the far-IR. 

In addition to a rising SED in the infrared due to this 
excess emission, Class I protostars display various ab- 
sorption features in their near- and mid-infrared spec- 
tra due to silicates and ices in the envelope. However, 
many parts of the infrared spectrum are not accessible 
from ground-based telescopes, which in addition often 
lack the sensitivi ty to observe faint objects. The Spitzer 
Space Telescope (|Werner et al.|[2004T ) opened a new win- 
dow onto low-mass star formation; with unprecedented 
sensitivity, the precursors of low-mass stars can be stud- 
ied in various evolutionary stages. 

Here we present 28 spectra from 5 to 36 fim 
of Class I objects in the Taurus star-forming re- 
gion observed wit h the Infrared Spectrograph 7 (IRS; 
iHouck et all I2004D on board of Spitzer. The 5-20 
fim spectra of 5 o f the 28 Clas s I ob jects were al- 
ready introduced in IWatson et alJ ([2004) ; these objects 
are IRAS 04016+2610, 04108+2803B, 04181+2654 B, 
04239+2436, and DG Tau B. Our IRS observations are 
part of a larger Infrared Spectrograph guaranteed-time 
observing program; we presented the Class II and III 
objects from our sample of 150 young stellar objects in 
Taurus in an earlier paper (jFurlan et al.ll2006h . 

Since the Taurus- Auriga star-forming r egion is young, 
(1-2 Myr) and relatively n earby (140 pc; iKenvon et al.l 
119941* feertout et afl ll999). many of its objects have 
been fairly well-studied at wavelengths covering almost 
the entire electromagnetic spectrum. Here we use SEDs 
constructed with data from the literature and our IRS 
spectra to constrain envelope models and to determine 
which parameters describe the circumstellar environment 
of the Class I objects in our sample. 

This paper is structured as follows: in § [2] we lay out 
our Spitzer IRS observations and data reduction; in § [3] 
we present the IRS spectra and SEDs of the 28 Class I 
objects in our sample; in § [4] we describe the envelope 
models used to fit 22 of our Class I objects and give 
detailed descriptions of the individual model fits; in § [5] 
we discuss the different appearances of Class I objects 
and give a summary of our models; and in §[6] we present 
our conclusions. 

2. OBSERVATIONS AND DATA REDUCTION 

Our Taurus targets were observed with the IRS on 
Spitzer during IRS observing campaigns 3 (2004 Febru- 
ary 6 to 8), 4 (2004 February 27 to March 5), 12 (2004 
August 30 to 31; only 3 targets), and 24 (2005 Septem- 
ber 12; only IRAS 04166+2706). All targets were ob- 
served over the full IRS range from 5 to 40 fim by using 
either the two low-resolution IRS modules (Short-Low 
[SL] and Long-Low [LL], 5.2-14 (im and 14-38 fim, re- 
spectively, A/AA ~ 90) or the SL module and the two 
high-resolution modules (Short-High [SH] and Long-High 
[LH], 10-19 fim and 19-37 fim, respectively, A/AA ~ 

7 The IRS was a collaborative venture between Cornell Univer- 
sity and Ball Aerospace Corporation, funded by NASA through 
the Jet Propulsion Laboratory and the Ames Research Center. 



600), depending on the object's expected mid- infrared 
flux. Thus, roughly half of our Class I objects were ob- 
served with SL and LL, and the other half with SL, SH, 
and LH. Our observations were carried out either in IRS 
staring mode, where the target is placed on two nod po- 
sitions along the slit, at 1/3 and 2/3 of the slit length, 
or in mapping mode, with 2x3 step maps on the target 
(3 steps in the dispersion direction and 2 steps in the 
spatial direction). More than half of the objects (19 out 
of 28) were observed in mapping mode. 

Data red uction was carried o ut using the SMART soft- 
ware tool (|Higdon et al.l 12004ft . Spectra were extracted 
after bad pixels, including so-called rogue pixels, were 
fixed by interpolation and after the sky in the low- 
resolution spec tra was subtracted using the off-order ob- 
servation (see IFurlan et al.l 120061 for more details) . A 
variable-width column extraction was performed for low- 
resolution spectra, while high-resolution spectra were ex- 
tracted using a full-slit extraction. The spectra were cali- 
brated by dividing them by a calibrator spectrum (a Lac 
(Al V) for low-resolution and £ Dra (K2 III) for high- 
resolution observations), and by mult iplying the result b y 
the calibrator's template spectrum ([Cohen et al.| [2003). 
Final spectra were obtained by averaging the two nod 
positions (for staring mode observations) or the two cen- 
tral map positions (for mapping mode observations). We 
also truncated all SH spectra below 13 fim and rebinned 
all high- resolution spectra to a resolution of A/AA = 200 
to achieve a more uniform representat ion of the spectra. 

As mentioned in lFurlan et alJ (|2006l ). we applied scalar 
corrections to each LH spectrum to account for the ab- 
sence of sky subtraction in the high-resolution mod- 
ules. Some minor adjustments (< 10%) were also nec- 
essary for SL to match it to the flux level of SH. 
For IRAS 04154+2823, 04158+2805, 04325+2402, and 
04368+2557, a multiplicative factor of 1.2 had to be ap- 
plied to SL, while DG Tau B is the only Class I object 
in our sample to require an even larger scaling: a multi- 
plication by 1.35 of SL2 (5.2-7.5 fim) and by 1.21 of SL1 
(7.4-14 fim). 

The flux mismatches between short- and long- 
wavelength modules can be attributed to the extended 
nature of our targets; even though all objects in our 
sample appeared as point sources when observed with 
the IRS, we know that most of our targets are some- 
what extended due to the presence of an envelope (e.g., 
iPark fc Kenvonll2002f) . There are two effects caused by 
the extension of the envelope influencing our IRS obser- 
vations. 

First, the LL and LH slits, which are 10.6" and 11.1" 
wide, respectively, cover a larger area and thus more of 
the envelope emission than the SL and SH slits, which 
are 3.6" and 4.7" wide, respectively. However, the cold, 
outer envelope regions do not contribute much flux at 
the IRS wavelengths; at a distance of 140 pc, the SL 
slit encompasses all emission within 250 AU from the 
central source, while LH includes emission out to 780 
AU, sufficient to sample almost all of the material con- 
tributing to the IRS spectrum. Only a few objects 
probably suffer from slit losses beyond about 30 fim, 
evidenced by a change in slope, and thus decrease in 
flux level, around 30 /im; this effect is most notable 
in IRAS 04108+2803 B, 04158+2805, 04181+2654 A, 
04239+2436, and 04295+2251. 
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Second, and more importantly, is the fact that proto- 
stellar envelopes are usually not spherically symmetric 
and that each of the IRS modules sample different re- 
gions of a protostellar system, since the slits of the var- 
ious modules are not aligned parallel to each other and 
are also of different sizes (see, e.g., Fig. Q}. For exam- 
ple, a protostar seen edge-on might show a dark central 
lane representing the disk and symmetric reflection neb- 
ulae representing inner envelope regions; if SL is aligned 
along the disk, LL will be oriented perpendicular to the 
dark lane and will sample more of the envelope emission. 

The larger scalar correction factors applied to some of 
Class I objects in our sample (especially in the case of DG 
Tau B) are likely due to this second effect. In particular, 
the centering and orientation of the narrow SL slit could 
result in some of the envelope emission being missed by 
the SL observation (see Fig. |T|) . However, in most cases 
only a small fraction of the warm dust contributing at 
IRS wavelengths is not included in the IRS observations, 
and therefore the spectra we obtained are believed to be 
a good representation of the mid-infrared emission of the 
inner envelope regions. 

Our absolute spectrophotometric accuracy is of the or- 
der of 5-10%, depending on the module used and the 
pointing uncertainty. The relative spectral accuracy 
is module-dependent; spectral features above the noise 
level are likely real in SL, LL2 (14-21 /urn) and SH, while 
in LL1 (20-36 fim) and especially in LH artifacts from 
unresolved calibration issues often do not allow reliable 
identifications of spectral features. Also the wavelength 
region where SL and LL meet (13.5-14.5 /im) is affected 
by this higher uncertainty due to possible artifacts at the 
order edges. 

Some of the Class I objects we observed with IRS are 
members of multiple systems (see Table [T]). Except for 
GV Tau A and B, they either are at separations that 
place all members well within the IRS slits, or their sep- 
arations are large enough to ensure that only one compo- 
nent is observed. Since the infrared emission from proto- 
stars is dominated by envelope emission, the multiplicity 
of the central source should have no impact on the mid- 
infrared spectrum, except when the configuration of the 
system is such that it influences the emission from the en- 
velope, e.g., by truncating it or causing a departure from 
spherical symmetry in the illuminating radiation field. 

We encountered reduction problems in a few objects; 
they are listed below. 

04108+2803B. This object is part of a mispointed 
mapping-mode observation; by using the spectra from 
the central map position for SL and SH and scaled spec- 
tra from the third map position for LH, we created a 
continuous IRS spectrum. This stitching adds a higher 
uncertainty to the absolute flux level in our spectra, but 
the overall shape of the spectrum should be only mini- 
mally affected. 

GV Tau (A,B). This binary is only separated by l'/3; it 
was observed with SL, SH, and LH in mapping mode, 
and it was not resolved in any of the modules. Since 
this object consists of two point sources with likely some 
extended emission surrounding them, their SH and LH 
calibration presented some problems. The final spectrum 
was composed by using those nod and map positions for 
each module in which the flux was highest, and in ad- 
dition SH was scaled to match up with SL. Due to this 



procedure, the absolute flux level of this spectrum car- 
ries a higher uncertainty. In fact, we multiplied the entire 
spectrum by 0.85 to match fluxes at 11.7 and 17.9 /im 
with our, yet unpublished, ground-based measurements 
at these wavelengths. 

HH 30. This object was too faint in SL2; the spectrum 
could only be extracted in SL1 and in LL. 
LkHa 358. The peak-up arrays contained bright to satu- 
rated sources during the SL observations, which resulted 
in increased noise in the SL part of the spectrum. The 
LL observation included 3 sources, LkHa 358, HL Tau, 
and XZ Tau, but since the observation was centered on 
LkHa 358, the other two sources entered only partially 
in the slit. In addition, all three sources were separated 
sufficiently in the spatial direction to allow extraction of 
the individual spectrum of LkHa 358. 

3. THE SAMPLE 

3.1. IRS Spectra 

Our Taurus sample is ba sed on objects analyzed by 
iKenvon &: Hartmannl (|1995l ) (from here on KH95); a 
large fraction of our Class I targets has identifiers from 
the IRAS Point Source Catalogue, since these objects 
were first identified as possible embedded protostars in 
the IRAS su rveys and confirmed i n follow-up obser- 
vations (e.g.. iBeichman et alj [l986t iMvers et all Il987t 
IKenvon et al.l Ii"990f ) . Some objects that are part of 
our Class I sample are identified as Class II objects in 
KH95: these objects are IRAS 04154+2823, 04158+2805, 
04278+2253, CoKu Tau/1, HL Tau, IC 2087 IR, and 
LkHa 358. We base our identification of these sources 
as Class I objects on the existence of extended, near- 
infrared nebulosity and possibly outflows around these 
objects, or the presence of a silicate absorption feature 
and a positive spectral index over the IRS range, or the 
detection of ice absorption features in the 5-8 /<m region 
as well as the CO2 ice feature at 15.2 /im. Therefore, this 
sample includes all Taurus objects we observed with the 
IRS that are thought to be at the protostellar evolution- 
ary stage and therefore surrounded, to varying extents, 
by envelopes, even though in a few cases the SED slope 
does not lie in the typical range for Class I objects. 

Figure [2] shows the IRS spectra of the 28 Class I ob- 
jects in our sample. Most have SEDs that are rising over 
the IRS spectral range, and almost all of them show nu- 
merous absorption features in their mid- infrared spectra, 
typical of embedded protostars (e.g.. lAlexander et al.1 
12003 iGibb etld][200l IWatson et al.ll2004h : wide bands 
due to amorphous silicates at 9.7 and 18 /tm (the lat- 
ter one is much less prominent), and features at 6.0 
/tin due to H2O ice, at 15.2 /im due to CO2 ice, and 
at 6.8 /tm due to organic compoun ds, possibly CH3OH , 
NHj~, or strongly po lar H2O ice ()van Dishoeckl 12004 
lAlexander et alj 120031 ) . In addition, the silicate absorp- 
tion feature at 9.7 /tm often displays a wider absorp- 
tion shoulder at the longer-wavelength end due to the 
broad H2O ice libration band at ~ 11-13 /tm. The long- 
wavelength wing of the CO2 ice feature also suggests that 
water ice is responsible, likely in t he form of a mix with 
C0 2 ice (|Ehrenfreund et al.lll999t) . 

Some protostars show an obvious self-absorbed silicate 
feature at 10 /im (e.g., IRAS 04154+2823, 04181+2654 
B; see Fig. [5Ji), where the silicate emission from in- 
ner disk regions is absorbed in cooler envelope re- 
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gions, causing a feature that is more or less flat, 
with some emission and som e absorption characteristics 
([Mitchell fc Robinsonlll98lD . The self-absorbed features 
probably indicate that the optical depth of the envelope 
along the line of sight is not as large as for the cases 
where a clear silicate absorption feature is seen. 

IRAS 04248+2612, 04264+2433, and CoKu Tau/1 
(Fig. [2^) have a silicate emission feature, but are sur- 
rounded by nebulosities and have an overall rising 
SED in the infrared. Therefore, they are likely more 
evolved Class I objects seen through low-density en- 
velopes, which have mostly dispersed and thus reveal 
the central star and its accretion disk. Based on the 
geometry of their near-infrared reflection nebulae, IRAS 
04248+2 612 and CoKu Tau /1 are likely seen close to 
edge-on ([Padgett et al.lll999fl . 

3.2. Spectral Energy Distributions 

Most protostars are found near or are embedded 
within dense mole cular cloud cores (jMvers et al.l Il987t 
lOnishi et al.1 [1998). Due to high extinction mostly due 
to their envelopes, protostars are usually faint or even 
undetectable in the optical. In the near-infrared and in 
the optical, the infalling envelope often appears as a re- 
flection nebula, seen in light scattered by dust grains in 
the inner parts of envelope cavities (e.g.. [Tarnura et all 
H99l IPadeett et al J 119991 IPark fc Kenvonll2002l) . In ad- 
dition, C lass I objects are often extended in the sub- mm 
and mm (|Chandler et al.lll998HMotte fc Andrell2001l ): at 
these wavelengths, the thermal emission from cold dust 
in the envelope is detected. 

To show the emission from protostellar systems at 
all wavelengths, in particular how the emitted power 
is distributed at different wavelengths, we generated 
spectral energy distributions with data from the liter- 
ature. In Figure EH we show the SEDs of the Class 
I objects in our sample, constr ucted with optical to 
mid-IR photometry mostly from iKenvon fc Hartmannl 
(1995), the 2MASS J, H, and K s fluxes, the IRAS 
12, 25, 60, and 100 am fluxes, IRAC 3.6, 4.5, 5.8 
and 8 .0 /tm fluxes (Hart mann et al.l l2005: Lu hman et alJ 
2006), sub -mm and mm fluxes, where available, and 
the IRS spectrum. The IRAS fluxes are mostly from 
the IRAS Faint Source Catalog; for CoKu Tau/1, GV 
Tau, HL T au, and L1551 IRS 5 w e adopted the IRAS 
fluxes from IWeaver fc Jones! (|1992ft . and for the 12 and 
25 /jm fluxes of DG Tau B and H H 30 the measure- 
ments from coadded IRAS scans bv ICohen fc Schwartd 
(|1987l ). The fluxes at 3 50, 450, 800, 850, 1100, a n d 130 
jim are mostly from iMoriartv-Schieven et alJ ([19941) . 
[Andrews fc WUliamsl (|2005f) . lYoung et all (|2003f) . and 
iMotte fc Andrei (|2001), while the few measurements at 



160 /im were taken bv lLadd et al.l ( 19911 ) with the Kuiper 
Airborne Observatory (see Table [3] for a list of the long- 
wavelength data used in the SED plots) . The fluxes were 
not corrected for reddening, since the high extinction of 
these objects is very likely intrinsic. 

The SED plots show that most of the power of a Class 
I object is emitted in the mid- to far-infrared; therefore, 
for many objects the IRS spectrum lies close to the peak 
of the emission. For most objects, the SED is either rising 
or flat over the IRS spectral range, with a deep silicate 
absorption feature at 10 /im as the most notable spectral 
feature. Only a few objects have SEDs that seem to peak 



at ~ 5 /im and decrease over the IRS spectral range; they 
might be more evolved protostars and will be discussed 
in § 15X51 

3.3. Notes on Individual Objects Excluded in the 
Modeling 

3.3.1. Class I Objects 

GV Tau. This object, also known as Haro 6-10, 
consists of two components separated by 1"3. Like T 
Tau, the GV Tau system consists of a T Tauri star 
(component A) and an embedded infrared compan- 
ion (component B). The binary is surrounded by a 
circumstellar disk and an almost edge-on circumbi- 
nary envelope, wit h outflows oriented perpendicular 
to the disk plane ([Menard et alJ Il993f ). GV Tau B 
is likely resp onsible for the outflows measured from 
this source ([Chandler et al.1 Il998t iHogerheiide et all 
fl998h . In the near-IR, the brightness ratio of the 
components varies considerably over the timescale of 
a few years (iLeinert fc Haasl 119891: iMenard et~aTI 119931 : 
iKoresko et al.l 1997). It is likely that the B component, 
which is more deeply embedded, varies more than the 
A component in brightness and generates the deep 10 
/zm silicate absorption feature, as well as the various ice 
features, observed in its mid-IR spectrum (see Figure 

m- . . 

HH 30. Optical and near -IR HST images ( [Burrows et al.l 
119961 : ICotera et al.l l200lh reveal that this object is seen 
edge-on; the central star illuminates a flared circum- 
stellar disk, which obscures the protostar and is seen 
in scattered light. In addition, a narrow jet is emitted 
perpendicular to the disk plane. Its appearance is 
reminiscent of that of DG Tau B, which is also seen 
edge-on and emits a jet. IWood et al.1 (|1998l ) modeled the 
emission of HH 30 and concluded that, in addition to a 
disk, only a tenuous envelope is required. The presence 
of a diffuse infalling envelope, as well as outflows, 
indicates that HH 30 is a protostar, though likely more 
evolved. 

The IRAS fluxes of HH 30, determined by 
ICohen fc Schwartz! (|1987h from coadded IRAS scans, 
are highly uncertain, since they are contaminated by 
emission from two nearby sources, HL Tau and XZ Tau, 
which are much brighter. Our IRS spectrum shows that 
the IRAS fluxes were overestimated by 1.5 to 2 orders of 
magnitude (see Figure [3]). HH 30 is among the faintest 
objects in our Taurus sample; its IRS spectrum suggests 
the presence of a deep silicate absorption feature, but 
its SED is only poorly determined. 

3.3.2. Class I/II Objects 



As indicated in § 13.21 a few objects in our sample 
have SEDs that are not typical of protostars. IRAS 
04278+2253 and IC 2087 IR have clearly decreasing 
SEDs from 5 to 40 /jm (see Fig. [5Ji) and would there- 
fore be identified as Class II objects based on the slope 
of their infrared SED. Also IRAS 04154+2823 and LkHa 
358 have SEDs that seem to peak around 5 /an and do 
not display prominent 10 /im silicate absorption features. 

However, IC 2087 IR (also known as Elias 18) has ice 
absorption features in the 5-8 /im spectral region, a 10 
/im silicate absorption feature, as well as a prominent 
CO2 ice absorption at 15.2 /im. It is likely surrounded 
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by a circumstellar disk, probably seen edge-on given it s 
high visual extinction Ay ~ 18 (|Shuping et al.l 1200 If ). 
iBowev. Adamson. fc Yatesl (|2003h measured the 10 fira 
spectrum of IC 2087 IR; they noted a strong silicate ab- 
sorption feature on a decreasing photosphere- like contin- 
uum. The absorption features we detect in the IRS spec- 
trum of IC 2087 IR, and the fact that it is surrounded 
by a reflection nebulosity (IC 2087), suggest that it is 
actually a protostar, surrounded by an accretion disk 
and some envelope material. Further evidence for its 
earlier evolutiona ry stage are a molecular outflow (e.g., 
iHever et al.lll987l ) and the presence of two Herbig-Haro 
knots, which are not within the outflow, but are ori- 
ented along the axis of symmetry of the reflection nebula 
(|Gomez et al.lll997t ). It has been suggested that, since 
IC 2087 IR has a rather weak outflow, it could be in 
transition be tween the embedd ed protostellar and the T 
Tauri phase (jTegler et al.lll995f ). 

In contrast to IC 2087 IR, IRAS 04278+2253 shows no 
silicate absorption feature, only a weak CO2 ice feature, 
and likely gas absorption features in the 5-8 /jm wave- 
length range, most prominently an H2O band centered 
at 5.7 pm. This featu re, similar to what is observed in 
some FU Ori objects (|Green et al.ll2006t ). is thought to 
arise in the upper layers of hot inner disk regions that 
are heated by accretion. IRAS 04278+2253 is a 6'.' 8 bi- 
nary, whose primary component is of G8 spectral type 
and has a relatively high mass accreti on rate (6.6 x 10~ 6 
m O vr~ 1 : IWhite fc Hillenbrand 2004 ). The seconda ry is 
of K7/M0 spectral type (|White fc Hillenbrandll200l and 
thus clearly fainter than the primary. Given this object's 
relatively large mass accretion rate, it is likely still a pro- 
tostar, but since no extended emission has been detected 
around it, there is likely only little envelope material left. 

LkHa 358 forms a small cluster of stars together 
with HL Tau and XZ Tau (a flat-spectrum source, see 
iFurlan et al.ll20 06): all three objects are found within a 
radius of less than 1'. All of them are likely embedded 
to some extent; while HL Tau has an obvious absorption 
feature and an SED typical of a Class I object, LkHa 
358 only has a weak absorption feature at 10 /im, and 
XZ Tau has silicate emission features at 10 and 20 /im 
and an overall flat SED in the mid-IR. Since HL Tau is 
the brightest source of the group, it is likely the dominat- 
ing contribution to the IRAS source 04287+1807; there- 
fore, the IRAS fluxes for LkHa 358 are not known. The 
IRS spectrum shows that LkHa 358 has a silicate feature 
that could be partly self-absorbed, as well as some weak 
ice absorption features. The late spectral t ype of LkHa 
358 ( M5.5 according to lKenvon et al.lfl998l and iLuhmanl 
20.00) suggests that this object will end up as a very low- 
mass star, with a mass close to that of a brown dwarf. 

IRAS 04154+2823 is very faint and could be a more 
evolved protostar; in fact, its 10 /im silicate feature seems 
to be self-absorbed, with a broad absorption centered at 
9.7 /im and emission between 10 and 12 /im, which is an 
indication of a tenuous envelope. Its SED can be fitted 
reasonably well with a low-density envelope model (see 
Figure [5]). 

4. ENVELOPE MODELS 

4.1. Model Description 

In this paper, we generated envelope models for 22 of 
the 28 Class I objects we observed in the Taurus star- 



forming region using the IRS on Spitzer to characterize 
some of the envelope parameters that result in the ob- 
served SED shapes. These models include the contribu- 
tions of the central protostar, the disk, and the envelope 
to compute the resulting SED of a protostellar system. 
We used a modeling code based on two types of mod- 
els that differ i n the initial configuration: so -called TSC 
models ( after iTerebev. Shu, fc Cassenl 11984 ). which as- 
sume an initially spherically sym metric molecular cloud 
core, and sheet-collapse models ([Hartmann et al.l [T994J : 
lHartmann. Calvet. fc Bossl Il996f ). which simulate the 
collapse of elongated cloud cores. Both models result 
in flattened density distributions, in which the emergent 
SED is strongly dependent on the inclination angle of the 
system. Below, we will briefly introduce these two types 
of models. 

4.1.1. Envelope Component 

In the TSC models, the radiative equilibrium tempera- 
ture is determined first by using the angle-averaged den- 
sity distribution for the infall region (p oc r~ 3 / 2 for 
r 3> R c , p oc r -1 / 2 for r -C R c ', R c , the centrifugal ra- 
dius, is where material falling in on the equatorial plane 
lands and thus represents the outer disk radius). In this 
region, departures from spherical symmetry become ap- 
parent only at small distances, when material falls onto 
a disk; at large distances, the TSC model can be de- 
scribed by spherical free-fall. After the temperature dis- 
tribution has been calculated, t he flattened, axially sym- 
metric density distribution from lTerebev. Shu, fc Cassenl 
(I1984T) is used to compute the flux emitte d by the sys- 
tem ([Kenvon. Calvet. fc Hartmannlll993al from here on 
KCH93). The models also include isotropic scattering. 
The source luminosity is not set according to the mass 
infall rate, which means that the disk accretion rate can 
be different from the mass infall rate of the envelope. 

Outflow cavities along the poles can also be included 
in t he models by specifying a cavity semi-op ening an- 
gle ([Kenvon et alj Il993bt (Calvet et all Il994f ). Inside 
the bipolar cavities, material has been evacuated; out- 
side, material will fall onto a disk with an inner radius 
Rmin = Rc sin 2 9 , where 9 is the half opening angle 
of the cavity, i.e. the angle between the polar axis and 
the streamline of an infalling particle at the edge of cav- 
ity at large distanc es, and R c is the centrifugal radius 
(|Calvet et al.| [l994). Thus, the shape of the cavities is 
determined by the streamlines of the infalling material. 

The inclusion of a cavity in the density distribution 
of the envelope causes departures from spherical symme- 
try in the temperature distribution. Scattered light from 
these cavities adds flux at near- infrared wavelengths, but 
the cavities also remove hot, inner dust. Therefore, for 
larger cavities, the net result is a decrease in near- and 
mid-infrared flux when the inclination angle is larger 
than the opening angle of the cavity. However, when 
the inclination of the system allows the observer to see 
inside the cavity and thus the star an d inner disk directly , 
the near-infrared emission increases ([Calvet et alJll994h . 

Sheet-collapse models start from the collapse of 
sheets that are i nitially i n hydrostatic equilibrium 
(lHartmann et all 119941 : lHartmann. Calvet. fc Bossl 
1996). This modification to the spherical infall scenario 
produces partially evacuated cavities, which, when 
aligned along the line of sight, reveal the inner parts of 
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the envelope and disk, thus increasing the object's opti- 
cal and near-infrared emission. Since the polar regions 
of the sheet collapse before the equatorial regions, polar 
cavities, which are re gions where low-densit y material is 
falling in, will result ( Hartmann et al.lll994h . With sheet 
collapse, the opening angle of the cavity is connected to 
the infall and not to the outflows; ev entually, outflows, 
but also wide-angle stellar winds (|Arce fc Goodman 
2002), might easily remove the low-density material and 
fully clear out the cavities, so in the end both effects, 
sheet collapse and mass loss from a protostar, could be 
responsible for the formation of cavities. 

4.1.2. Disk Component 

The disk is modeled as an optically thick, steady ac- 
cretion disk irradiated as a flat disk by the central star 
(T(r) oc r 3 / ; T = T accretion + T irradiation ) , oriented 
perpendicular to the rotation axis of the envelope. As 
the outer disk regions, which emit mostly in the sub-mm 
and mm, are optically thin, and only the inner disk re- 
gions are optically thick, the outer disk radius was set at 
10 AU to represent the optically thick part of the disk. 
However, the disk is expected to extend out to i? c ; emis- 
sion from the protostellar envelope will be obscured if it 
lies behind the disk along the line of sight. Therefore, to 
take this obscuration effect into account, a second set of 
models was calculated with the disk radius set equal to 
R c . The latter set was used for the envelope and stellar 
emission, while the set with the smaller disk radius was 
used for the disk emission. In this way the disk com- 
ponent in our models represents only the optically thick 
disk regions; the optically thin outer regions would likely 
add small amounts of additional flux at sub-mm and mm 
wavelengths. 

This simple disk model also does not include an inner 
disk wall at the dust truncation radius, which would add 
flux at near-IR wavelengths, in particular for systems 
with higher inclination angles i, where the disk emission 
is reduced due to its cos(i) dependence. In addition, 
since the disk model does not include an optically thin 
disk atmosphere that generates a 10 jim si licate emis- 
sion feature (see, e.g.. iD'Alessio et all Il999f ). the shape 
of the silicate feature is only determined by the envelope. 
The material in the envelope will extinguish the disk and 
thus generate a 10 /zm absorption feature in the contin- 
uum emission arising from the disk, while the silicates 
in the envelope will produce an emission feature when 
the density is low and an absorption feature for higher 
densities. Thus, the silicate feature might not be mod- 
eled realistically for Class I objects with low envelope 
densities, but should be a good approximation for more 
typical, higher-density envelopes. 

The inner disk radius is set by the star's magneto- 
sphere, which channels the accretion flow onto the star, 
and is roughly at the corotat ion radius, at a few stellar 
radii from the stellar surface (|Shu et al.lll994l ). However, 
in some models the inner disk radius was set equal to the 
stellar radius, implying a boundary layer. The accretion 
luminosity is set as the sum of the disk accretion luminos- 
ity, generated as matter accretes through the disk, and 
of the hot component luminosity, produced by the dissi- 
pation of the accretion shock on the stellar surface. For 
the cases where the inner disk radius is equal to the stel- 
lar radius, half of the accretion luminosity is emitted by 



the disk, and the other half by the boundary layer. For 
larger inner disk radii, the accretion shock accounts for 
a larger fraction of the accretion luminosity than matter 
accreting through the disk. 

4.1.3. Dust and Ices 

The dust in the envelo pe is composed of s mall sili- 
cate and graphite gr ains (|Draine fc Led Il984l ). troilite 
(Fe S), and water ice (IPollack et all 11994 ). as well C0 2 
ice ()Ehrenfreund et al.lll996f ) . The particle size distribu- 
tion follows a power law with particle radius a: n(a)da cx 
a~ 3 ' 5 da, with maximum sizes of 0.3 /im. This maximum 
grain size corresponds to t he typical value found in the 
interstellar medium (e.g., Q vlathis, Rumpl. fc Nordsieckl 
Il977t iWeingartner fc Drainel l200lh and is motivated by 
the notion that envelopes likely contain pristine material. 
Our data confirm this idea, since the shape of the silicate 
absorption feature in most of our Class I objects indi- 
cates that the dust is dominated by small grains (see also 
lOsorio e~a l. 2003). The abundance of the various com- 
po nents of the dust mix ture is similar to the one adopted 
by lOsorio et al.1 (12003 ) who adju sted the abundances 
given in iDraine fc Led (fl98l and IPollack et all ([1994) 
to fit a model to the spectrum of L1551 IRS 5. We do 
not include organics in our mixture; we assume that the 
carbon is present in the form of graphite only. The frac- 
tional abundances with respect to the mass of the gas for 
silicates, graphite, troilite, and H2O ice are ( s u = 0.004, 
Qra = 0.0025, Ctro = 0.000768, and (ice = 0.0005, re- 
spectively. To get a better fit to the observed SEDs, 
these fractional mass abundances could be varied until a 
good match was found, but in order to limit the number 
of free parameters in the fit, we kept the abundances of 
these components constant. Only the CO2 ice abundance 
was adjusted for each object to reproduce the depth of 
the CO2 ice feature. 

The opacities used in the current model should be re- 
garded as approximations: the water ice opacities used 
in the model do not generate a feature at 6.0 /zm, which 
is detected in our IRS spectra of Class I objects, and 
other ice opacities, e.g. for methanol ice, which might be 
responsible for the feature at 6.8 /im, are not included. 
Also, the 15.2 /im ice feature is treated as pure CO2 ice, 
but in reality is likely mixed with water and other ices , 
making the feature broader (jEhrenfreund et all fl999). 
Therefore, our envelope models are not able to reproduce 
the majority of the ice features, and the most prominent 
ice feature, the CO2 feature, can only be approximated 
in its depth, not its width. On the other hand, since the 
purpose of this study is a characterization of the over- 
all SEDs of Class I objects, and thus of their envelope 
structure, we will neglect the fact that the ice features 
will not be reproduced well or at all. 

4.1.4. Model Parameters 

The input parameters for the models are the total sys- 
tem luminosity {L sta r + L acc ), the fraction of the total 
luminosity that arises from the star (r] s t a r), the radius of 
the star (assumed to be 2 R©), the minimum disk radius, 
the centrifugal radius (which is essentially the outer disk 
radius), the inner and outer envelope radii, the flatten- 
ing of the envelope (for sheet-collapse models), and the 
reference density p\, which is the density the envelope 
would have at a radius of 1 AU for the limit R c = 0. The 
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latter parameter, p\, is rel ated to the mass infall rate Mj 
and th e central mass M* ([Kenvon. Calvet. &: Hartmannl 
Il993al) : 

Kq m-14/ j^j ^ f M * V 1/2 -3 

pi = 5.3 x 10 - — K~~r~r r — 71— gem 

^lO-SMgyr-iy V 1M o/ 

For a fixed central mass, pi determines the mass infall 
rate. The inner radius of the envelope is determined 
by the dust destruction, which, for typical dust made 
of silicates, is at the dust sublimation temperature of ~ 
1500 K. If a cavity is included, the inner envelope radius 
will be farther out. The outer envelope radius is con- 
strained particularly by sub-mm and mm observations, 
since this parameter affects the SED almost exclusively 
beyond about 100 pm; typical values are several 1000 
AU. 

For sheet-collapse models, the flattening of the enve- 
lope is described by r) = R ou t/H, where R ou t is the outer 
radius of the envelope and H is the scale height of the 
original sheet, an isothermal, infinite, flat layer in hy- 
drostatic equilibrium, and is defined b y H = cj/jnGT,), 
with S equal to the surface density (|Hartmann et al.l 
Hflfll . Over the course of the collapse of a flattened 
molecular cloud core, it is thought that rj increases, re- 
flecting an increase in the size of the cavity. 

The emergent SED is computed with the disk radius 
(set equal to 10 AU for the disk component, and to R c for 
the stellar and envelope components), the inclination an- 
gle to the line of sight, and the opening angle of the cavity 
as parameters. Scattering is also included, which causes 
an increase in flux at near- infrared wavelengths (< 3 pm). 
The parameters that affect a modeled SED the most are 
pi, R c , and the inclination angle. The higher the density, 
the deeper the silicate absorption feature and the longer 
the wavelength at which the SED peaks. A high incli- 
nation angle will also increase the silicate feature depth, 
in addition to decreasing the emission below 8 pm and 
increasing the emission in the far-infrared. The centrifu- 
gal radius, besides representing the outer disk radius, 
also defines the region where infalling material piles up; 
thus, for smaller R c high-density material can be found 
closer to the central protostar, resulting in more mid- 
to far-infrared emission and a deeper silicate absorption 
feature. A sketch of our model envelope, with important 
model parameters identified, is shown in Figu re [U 

From previous studies (e.g. jKenvon et aill993bf ). typi- 
cal values for young protostars in the Taurus star-forming 
region are an R c of ~ 100 AU, and mass infall rates a few 
times 10 -6 M Q yr _1 . Observations exclude centrifugal 
radii larger than about 1000 AU and mass infall rates 
larger than 10 -6 M Q yr^ 1 . From equation [1] it follows 
that typical val ues for p\ are of the order of 10 ~ 14 to 
10" 13 g cm- 3 . lHartmann. Calvet. h Bosi (|1996l ) find 
that characteristic values of 77 are between 2.5 and 3.5 for 
the main infall phase of the envelope, resulting in large 
opening angles for the cavities. Over time, r\ should in- 
crease, up to a maximum value of ~ 4, and the cavities 
become larger. Inclinations are usually larger than 60°, 
but this could be a selection effect, since the SEDs of 
low-inclination protostars with large cavities would re- 
semble the SEDs of T Taur i stars and therefore not be 
classified as Class I objects (|Kenvon et al.lll993b[ ). 

In the following section, we will introduce our envelope 



models for 22 Class I objects of our sample. In Figures[5] 
to l28l (excluding Figures l6l and fl"8|) . we show the SED of 
each object together with the best-fit model (determined 
by eye) and its components, the star (dotted gray line), 
the disk (dash- dotted gray line), both extinguished by the 
envelope, and the envelope (long-dashed gray line); the 
sum of all components is represented by the thick, solid, 
gray line. All model protostars were placed at a distance 
of 140 pc. 

4.2. Models of Individual Class I Objects 
4.2.1. 04016+2610 

iPadgett et al.1 (|1999f ) imaged this object in the near- 
IR with HST; it is surrounded by a nebulosity on one 
side, with the protostar at its base. This might indi- 
cate that we are observing IRAS 04016+2610 through 
the opening of a large cavity in the envelope. A bipolar 
molecular outflow is oriented along the axis of symmetry 
of the nebula (e.g.. lGomez et al]ll997tlHogerheijde et al.1 
19981): its inclinat i on an gle was estimated to be 60° by 
Hog erheiide et alj (jl998l ). 

In Figure [5] we show an envelope model for IRAS 
04016+2610 whose parameters were adjusted to yield a 
good fit of the near-IR to mm SED. A sheet-collapse 
model with r\ — 1 had to be adopted; TSC models did 
not yield enough flux on the short- and long- wavelength 
sides of the 10 pm silicate feature. The IRS spectrum is 
fitted reasonably well between about 8.5 and 25 pm; the 
discrepancies at the shorter wavelength side might be re- 
solved by adjusting the abundances of the dust and ices 
in the model, while the lower IRS fluxes beyond 25 pm 
could be an aperture effect, with the IRS slits truncating 
the regions contributing to the 20-40 pm region. 

The model shown in Figure [5l with a density of 
4.5 x 10~ 14 g cm -3 , an R c value of 100 AU, and an 
inclination angle of 40°, is one of several that yielded a 
comparable good fit. In particular, values for pi down to 
3.5 x 10~ 14 g cm -3 and for R c down to 70 AU, as well as 
inclination angles within 30°-50°, would yield acceptable 
fits. The lower values for p\ and R c would match previ- 
ous modeling results of KCH93, who derived 3.2 xl0~ 14 
g cm -3 and 70 AU, respectively, for these two parame- 
ters. They would also be r oughly consistent with results 
from lWhitnev et alJ <| 19971 ) . who derived an even smaller 
R c of 50 AU and an inclination angle of 57°. On the other 
ha nd, our best- fi t mod el agree s very well with the res ults 
of lEisner et all (|2005h and of iGramaio et al.l (j2007h . A 
much smaller R c would result in an SED that is too nar- 
row and contain too much flux in the 20-50 pm region, 
while a large R c would require larger inclination angles 
to match the 10 pm feature and longer-wavelength IRS 
flux, but it would also decrease the emission below 8 pm. 
IStark et all (|200l require a larger centrifugal radius of 
300 AU, combined with a larger value for p\ and an in- 
clination angle of 65°, to reproduce the near-IR images 
of IRAS 04016+2610. 

Our derived range for the inclination angle, 30°-50°, is 
somewhat smaller than the value of 60° derived from ob- 
servations (Hog erheiide et alj|1998t IPadgett et ai] [l999) 
and modeling (KCH93). In addition, our derived cav- 
ity semi-opening angle of 5° is smaller than ex pected 
from the images; for example, IStark et alj (|2006f ) derive 
a value of 25° for 0. If we adopted a larger cavity, the 
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emission from near-IR wavelengths to about 20 pm would 
decrease and therefore not fit the shape of most of the 
IRS spectrum. Finally, our adopted luminosity of 4.5 
Lq is larger than the bolometric luminosity of 3.7 L Q 
(KH95), but necessary to yield enough flux. 

4.2.2. 04108+2803 B 

This protostar is part of a binary; it is separated by 
about 21" from component A, which is less embedded, 
displays a silicate emission feature, but is about a factor 
of 5 fainter in the mid-infrared (see Figure [S]) . There- 
fore, component B clearly dominates in the mid-IR; while 
04108+2803 A probably is a reddened T Tauri star, 
04108+2803 B is a Class I object. Both source s have 
no optical cou nterparts (e.g . , iTamura et al .1 1 1 99lh . 

According to lYoung et all (j2003j), IRAS 04108+2803 B 
is a Class I object based on its bolometric temperature of 
179 K, but would be classified as a Class object based 
on its submillimctcr luminosity. This IRAS source is not 
resolv ed in the 450 and 850 pm maps of lYoung et alJ 
(l2003f): since it is als o very faint at 1.3 mm (39 mJy; 
iMotte fc Andrei [2001h and does not have an outflow, it 
is likely an evolved protostar with a smaller envelope. 

The envelope model for IRAS 04108+2803 B requires 
relatively low density (p\ = 1.5 x 10 -14 g cm -3 ), a small 
centrifugal radius (R c = 40 AU), a moderate-size cavity 
(0 = 10°), and an inclination angle of 40° (Figure [7|). 
The outer envelope radius was set at 6000 AU, which 
is larger than the projected binary separation of 3000 
AU, but necessary to fit the sub-mm and mm emission. 
This suggests that component A is likely obscured by this 
envelope; alternatively, IRAS 04108+2803 B could be 
surrounded by a smaller envelope, but a fairly large cir- 
cumstellar disk, which would add to the long-wavelength 
emission. 

The IRS spectrum is fitted rather well with this en- 
velope model; the change in slope of the spectrum be- 
yond 25 /im could be an effect of slit truncation. The fit 
of the 10 /im feature and of the slope of the IRS spec- 
trum between 10 and 25 /im suggests that the values 
for the density and the centrifugal radius are relatively 
we ll-constrained. Previo us models by KCH93 and also 
by IWhitnev et al.l (|1997f) obtained higher values for p\ 
and R c (about dou ble our values ) , but with comparable 
inclination angles. lEisner et al.l (|2005l ) derived a large 
value for pi, too, but an R c similar to our value; their 
best-fit model included a massive disk surrounded by a 
500 AU envelope. 

In our models, we can exclude R c values larger than 
about 50 AU, since they would result in a decrease in flux 
from 15 to 40 /im due to an overall shallower SED shape. 
A reference density up to 2.0 x 10~ 14 g cm -3 would still 
fit the observed SED, if combined with a smaller incli- 
nation angle (~ 25°) to reduce the depth of the silicate 
absorption feature. The cavity semi-opening angle of 10° 
is required to fit the IRS spectrum; a smaller opening an- 
gle would result in an increase in flux in the 10-25 /im 
wavelength range. The inner disk radius was set at 3 
stellar radii to decrease the flux in the 2-8 pm region; 
this parameter does not affect the remaining parts of the 
SED. 

4.2.3. 04154 +2823 



IRAS 04154+2823 might be a transition object be- 
tween the Class I and II stage, since it has an SED 
that is somewhat decreasing in the mid-IR, and it also 
does not display a prominent silicate absorption feature. 
However, its spectrum shows a clear CO2 ice absorp- 
tion feature at 15.2 /im, and its overall infrared SED 
is relatively flat, suggesting that some envelope mate- 
rial should be present . It is also very faint in the optical 
(jStrom fc Strom 1994), possibly indicating larger extinc- 
tion along the line of sight. It is somewhat reminiscent 
of IRAS 04158+2805, which is also faint and has a rela- 
tively flat SED over the IRS spectral range. 

This protostar has a low luminosity {Lboi = 0.3-0.4 
L Q ; see Table 01; its envelope model requires a very low 
reference density of p\ = 7.0 x 10~ 15 g cm~ 3 and a very 
small centrifugal radius of 10 AU to reproduce the mid- 
IR spectrum as well as the sub-mm measurements (Fig- 
ure [8]). However, this model overestimates the flux be- 
tween 20 and 100 pm and slightly underestimates the 
mm flux. A luminosity value higher than about 0.4 L© 
would result in an increase in the long-wavelength flux, 
which would fit the mm data better, but yield even more 
flux between 20 and 60 pm.. An increase in p\ would 
have a similar effect; in addition, values of p\ larger than 
about 8.0 x 10 -15 g cm' 3 can be excluded given the con- 
straint of the weak silicate absorption feature. The low 
inclination angle (i = 20°) is also well-constrained by the 
10 pm feature and the shape of the SED. 

4.2.4. 04158+2805 

According to M enard et alJ (|2001h who imaged this 
source in the optical, IRAS 04158+2805 is likely sur- 
rounded by an edge-on disk and a bipolar nebula; a jet 
is oriented perpendicul ar to the disk p l ane, i n a geom- 
etry similar to HH 30. iPark fc Kenvonl (|2002t ) observed 
that this object is not extended in the near-IR or at mil- 
limeter continuum wavelengths, and therefore it could 
be a heavily re ddened T Tauri star and not a Class I ob- 
ject. Recently, lAndrews fc Williams! (|2007f) modeled the 
SED of IRAS 041 58+2805 with only a large disk com- 
ponent. However, iKenvon et alJ (|1998h concluded from 
their optical spectra and optical and near-IR photometry 
that IRAS 04158+2805 is a Class I source. Unpublished 
HST images at optical and near-IR wavelengths (pro- 
grams 9103 and 10603, PI K. Stapelfeldt and D. Padgett, 
respectively) show an extended, one-sided conical nebu- 
losity, confirming that this source is surrounded by an 
envelope. We note that th is object is of a late spectral 
ty pe (M5-M6 according to lLuhmanl l2006; M6 according 
to I White fc Hillenbrandll2004fk it will probably become a 
very low-mass star, close to the hydrogen-burning mass 
limit. 

The SED of IRAS 04158+2805 can be fitted with 
an envelope model with relatively low density (p\ = 
2.0 x 10~ 14 g cm -3 ), small centrifugal radius (R c = 60 
AU), small cavity (6 = 5°) and an inclination angle of 
30° (Figure [9]). The inner disk radius was set at 3 stellar 
radii to decrease the amount of flux emitted at near-IR 
and mid-IR wavelengths out to 8 pm. The outer enve- 
lope radius was set at 5000 AU to somewhat decrease 
the emission in the sub-mm and mm. Nevertheless, the 
model overestimates the flux at wavelengths longer than 
30 pm, but lies in between the two measurements around 
1 mm. The discrepancy between the two long- wavelength 
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data points could be explained by the large difference of 
the two beam sizes (see Table [2]); the 880 /jm measure- 
ment likely misses some extended emission. 

We adopted a luminosity of 0.3 L , which is 50% larger 
than the bolometric luminosity we measured by integrat- 
ing the fluxes fluxes of the SED; a model with L = 0.2 
Lq, a smaller density (~ f.O x 10~ 14 g cm -3 ), and a 
centrifugal radius around fO AU would yield a compa- 
rable fit, but still overestimate the emission starting at 
20 /an. Overall, the model parameters are mainly con- 
strained by the IRS spectrum; the depth of the silicate 
feature and the shape of the IRS spectrum from 12 to 
30 pm constrain the reference density p\ and the cavity 
semi-opening angle, respectively, reasonably well. With 
the current reference density, a much smaller centrifugal 
radius would result in even higher fluxes in the 30-60 
pm region, while a value larger than about 90 AU would 
decrease emission from about 15 to 40 pm. A pi value 
of about 1.5 xlO -14 g cm -3 , combined with an R c value 
of 50 AU and a slightly larger inclination angle of 40° 
would also yield an acceptable fit. We can exclude an 
edge-on orientation based on the shape of the IRS spec- 
trum. Only the availability of accurate, far-IR to mm 
fluxes will improve the model fit of IRAS 04158+2805. 

4.2.5. 04166+2706 

This IRAS source is very faint in the near-IR : no so urce 
is detected in 2M ASS but both lKenvon et all (j 19901 ) and 
iPark fc Keiivonl (|2002D detected a faint near-IR source 
close to the IRAS position, which we adopted for the 
pointing of our IRS observations. The latter authors de- 
scribe the appearance of I RAS 04166+2706 as low sur- 
face brightness neb ulosity. iMotte fc Andrei (|2001[ ) and 
Tafall aet al.l (|2004l ) detected a dense, round core at mm 
wavelengths coinciding with the IRAS source. Thus, this 
Class I object is likely deeply embedd ed. It also has bipo- 
lar outflows (jBontemps et al.1 |l996) , whi ch are highly 
collim ated and of extremely high velocity (tTafalla et al.l 
2004), consistent with this object's early evolutionary 
stage, possibly placing it at the Class stage. The 
high degree of sym metry of the red- a nd blueshifted 
gas of the outflows (Tafa lla et alj I2004D suggests that 
this object is likely seen close to edge-on. For the SED 
plot, we included th e H and K magnitudes measured by 
iKenvon et all (|199Ct ). 

The SED shape of IRAS 04166+2706 further indicates 
that it is likely seen at a large inclination angle, close 
to edge-on. In fact, an envelope model that fits it well 
(see Figure ITUl) has an inclination angle of 85° . Smaller 
inclination angles (70°-80°) would result in models that 
match the observed spectrum from 15 to 30 /jm more 
closely, but they would overestimate the emission from 
the near-IR to 8 fim. Increasing the cavity semi-opening 
angle by just 1° from its current value of 6° would lower 
the flux between 8 and 30 pm, resulting in a better fit of 
the slope of the spectrum beyond 12 pm, but a poorer fits 
of the silicate feature. An inner disk radius of 5 stellar 
radii was adopted to lower the flux below 8 pm. 

A reference density of 4.5 x 10~ 14 g cm -3 and a large 
centrifugal radius (R c — 300 AU) are required to fit the 
depth of the silicate feature and also the width of the 
SED. Values of p x and R c down to 3.5 x 10" 14 g cm" 3 
and 200 AU, respectively, would also yield good fits; only 
the silicate feature would be somewhat narrower. The 



best-fitting model of KCH93, who had fewer constraints 
available, resulted in p\ = 1.0 x 10~ 13 g cm -3 , R c =70 
AU, and i — 30°; a high- inclination model with the same 
R c , but i = 90° and pi = 3.2 x 10~ 14 g cm" 3 was also 
suggested as a fit. This underlines the importance of the 
IRS spectrum to discriminate between different models. 

We note that our model does not fit the 850 and 1300 
pm measurements, which were done with a larger beam 
size (see Table [2]), but it fits the data taken at similar 
wavelengths, but with ~ 20" apertures. This could indi- 
cate that some extended, cold dust outside of the infall 
regi on contributes to the long - wavelength emission (see 
also [Tayawardhana et al. 2001]). 



4.2.6. 04169+2702 

This object, which is not detected at optical wave- 
lengths, is surrounded by a cometary reflection nebula 
in the near -IR and associated with Herbig-Haro e mis- 
sion knots (jTamura et "all 119911; iGomez et alJ 119971) . It 
also has bipolar outflows (e .g., iMoriartv-Schieven et al.l 
1992; Bontemps et al.l[l996| ), which are perpendicular to 
an elongated envelope structure, about 2200 AU x 1100 
AU in size and incl ined at 60° from the line of sight 
(|Ohashi et al.lll997rj ). Measurements of continuum emis- 
sion at 2.7 mm did not resolve the source on scales of 
2"-3", probably indicating that this emission ari ses from 
a disk of radius < 150 AU (|Ohashi et al.lll997bl ). 

The SED of IRAS 04169+2702 is fitted well by a sheet- 
collapse envelope model (r) = 1) with moderate density 
(pi = 3.2 x 10~ 14 g cm -3 ), a centrifugal radius of 100 AU, 
basically no cavity and a high inclination angle (Figure 
ITT]) . Compared to the results obtained by KCH93, both 
our reference density and centrifugal radii are smaller 
by about a facto r 3, and our inclinati on angle is larger 
(75° versus 30°). I Whitney et~aH ([19971 ) modeled this ob- 
ject with an R c of just 10 AU, a reference density half 
our value, a cavity semi-opening angle of 11°, and an 
inclination angle of 45°. However, our model is better 
constrained. Both r] and the cavity semi-opening angle 
are determined by the shape of the IRS spectrum be- 
yond about 12 pm. In particular, a TSC model would 
yield less flux from about 10 to 100 pm given the same 
envelope parameters. The large R c value is required to 
fit the long-wavelength part of the IRS spectrum, while 
the value for the reference density is constrained by the 
shape of the silicate absorption feature and the peak of 
SED. The 10 pm feature and the SED shape also con- 
strain the high inclination angle, which has to be larger 
than 70° to yield a good fit. 

The luminosity used for the model (L = 1.5 L Q ) is 
twice as high as the bolometric luminosity of 0.8 L mea- 
sured by KH95, but close to the value we calculated by 
integrating the fluxes of the SED, 1.2 L^, and t o the 
value of 1.4 L Q determined bv lKenvon et al.l (|1990f ). We 
can exclude a lower luminosity from our model fit, since 
it would not yield enough emission to match the observa- 
tions. As for IRAS 04166+2706, our model does not fit 
the long-wavelength emission measured with large aper- 
ture sizes, but it matches the observations taken with 
smaller apertures. Again, we suggest that some cold dust 
just outside the infalling envelope is responsible for the 
extended mm emission. 

4.2.7. 04181+2654 A 
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This IRAS source is a member of a wide binary sys- 
tem; component B lies to the northwest of component 
A, at a projected distance of 31", which corresponds to 
about 4300 AU at the distance of Taurus. In addition, 
the IRAS source 04181+2655 lies somewhat less than 
1' to the northwest of 04181+2654 B. A bipolar out- 
flow is centered on IRAS 04181+2655, but no outflow 
has been detected from the IRAS sources to the south 
(jBontemps et al.l [199 61. which could be the result of a 
low inclination angle along the line of sight of these two 
objects. 

Of the two components of the 04181+2654 binary, the 
B component seems to be more deeply embedded, since 
it is fainter than A by a factor of 12 in the J band, 
but only by about a factor of 2 in the K band. Over 
the mid-IR spectral range, component B is fainter than 
component A by a factor of 2 to 3. The silicate feature 
of IRAS 04181+2654 B is seen partly in self-absorption, 
indicating that it is oriented along the line of sight such 
that part of the silicate emission originating in the inner 
d isk is absorbed in the en velope. 

iMotte fc Andre! (|200lD measured the 1.3 mm flux 
of IRAS 04181+2654 A, which is dominated by emis- 
sion from an exten ded envelope. According to 
iPark fc Keiivonl (|2002D . both A and B are extended at 
mm continuum wavelengths, but not in the near-IR. 
They suggest that both sources are protostars; there- 
fore, they are likely surrounded by an envelope which 
could be tr uncated in its outer parts. Images at 0.9 /tm 
obtained bv lEisner et alJ (|2005h reveal a faint nebulosity 
around IRAS 04181+2654 A, possibly tracing scattered 
light from an envelope cavity; the B component is not 
detected. 

The SED of IRAS 04181+2654 A can be fitted by an 
envelope model with a reference density of 2.0 x 10~ 14 
g cm~ 3 , a centrifugal radius of 50 AU, a cavity semi- 
opening angle of 7°, and an inclination angle of 20° (Fig- 
ure [T^]). In order to fit the SED shape from the near-IR 
to 8 /im, we set the inner disk radius to 5 stellar radii; 
even with this value, the near-IR fluxes are somewhat 
underestimated. We chose an outer envelope radius of 
10000 AU; such a large envelope would encompass com- 
ponent B, too, but is necessary to yield enough flux in 
the mm wavelength range. The fact that the model still 
underpredicts the flux at 1.3 mm suggests that the mm 
measurement probably includes some contribution from 
an envelope around component B. 

The density is constrained by the depth of the sili- 
cate absorption feature, which is rather shallow; it could 
be slightly larger, but then require a somewhat smaller 
cavity semi-opening angle to decrease emission in the 8- 
25 yum range. The slope of the spectrum from 12 to 
30 /im constrains both the opening angle of the cav- 
ity and the value of R c . A larger R c would result in 
a broader SED, but at the same time decrease the flux 
levels in the long-wavelength part of the IRS spectrum. 
A smaller R c would generate a steep rise of the flux levels 
beyond 12 /mi. KCH93 modeled both IRAS 04181+2654 
A and B with an envelope with a higher reference den- 
sity, 3.2 x 10~ 14 g cm~ 3 , a somewhat larger R c (70 AU 
versus our 50 AU), a nd a comparable inclination angle. 
IWhitnev et al.l (|1997h used the same value for R c as we 
did in our model, but with almost double the values for 
reference density and inclination angle. 



Overall, the model for IRAS 04181+2654 A is poorly 
constrained, given the lack of sub-mm data points and 
the uncertainty of the 1.3 mm measurement due to the 
presence of component B. We did not model that com- 
ponent, since there are no sub-mm or mm measurements 
available that would constrain the long-wavelength part 
of the SED. 

4.2.8. 04239+2436 

In the near-infrared, a symmetric reflection nebula sur- 
round s this IRAS source, which is a binary separated by 
0'.'3 (jReipurth et al.ll2000h . The brighter component (at 
2.2 /im) of this sys tem drives a large, asy mmetric Herbig- 
Haro jet, HH 300 (|Reipurth et al.ll2000f) : the asymmetry 
suggest that the inclination angle of this object is closer 
to face-on than edge-on. Maps of this Class I object 
at 450, 850, 1100, and 1300 /im show that it is a point 
source with some diffuse emission from circu mstellar ma- 
terial ([ChiEEiraLl HiHIl msnereEaDEOOl). Therefore, 
IRAS 04239+2436 is exte nded from infrared to millime- 
ter wavelengths (see also lPark fc Kenvonll2002T ). typical 
for Class I objects. 

The SED of IRAS 04239+2436 requires a model with 
relatively low density (pi — 1.3 x 10~ 14 g cm -3 ) and 
very small centrifugal radius (R c = 10 AU) (Figure [13]). 
A higher density would increase the depth of the silicate 
feature, while a larger centrifugal radius would decrease 
the emission in the 10 to 40 /im region and result in a 
shallow silicate absorption feature. Increasing p\ and R c 
up to about 2.0 x 10 -14 g cm~ 3 and 20 AU, respectively, 
would still result in good fits if the luminosity was in- 
creased by about 20%. 

Previous models by KCH93 yielded a larger density, 
3.2 xl0~ 14 g c m" 3 , and a centrifuga l radius of 70 AU, 
while models by IWhitnev et al.l (119971) resulted in values 
for p\ and R c very similar to ours. Eisner et all (|2005l ) 
derived a reference density almost twice our value and 
a centrifugal r a dius of 30 AU; more recent results by 
iGramaio et alJ (120071) found the same centrifugal radius 
as lEisner et alJ (|2005l ). but a reference density similar to 
ours. Our inclination angle of 15° is smaller than the 
values derived by previous modeling efforts, which range 
from 30° to 63°; a value of up to i = 30° would yield an 
acceptable fit, albeit with somewhat decreased emission 
from the near-IR to about 25 /tm. The low inclination 
angle is well-constrained by the SED, and it is also con- 
sistent with the near-IR morphology and observed asym- 
metry of the Herbig-Haro jet (see, e.g., Rcipurt h et al.l 
f2000f) . 

With the current values for R c and p\ , the model does 
not reproduce the shape of the silicate absorption fea- 
ture very well, especially at the long-wavelength side; 
this could be an effect of dust and ice composition, since 
the overall width of the SED and shape of the IRS spec- 
trum beyond 15 /im suggest such a small R c value. 

4.2.9. 04248+2612 

This IRAS source is actually a subarcsecond binary 
whose components are roughly similar in brightness. 
Its IRS spectrum shows a silicate emission feature 
and some very weak ice absorption features. In the 
near-IR, it is surrounded by a bipolar reflection neb- 
ula with a somewhat complex morphology, but likely 
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seen close to edge-on (Pa dgett et al.l [l999h . This ob- 
ject is also somewhat extended at 450 and 850 ^m; 
the extension parallel to the bipolar nebulosity seen in 
the near-IR is likely tr ac ing dust heated in an out - 
flow (|Young et al.l [2003h . IWhite k Hillenbrand! IpOOl ) 
note that IRAS 04248+2612 is the lowest-mass protostar 
known that drives a molecular outflow and an HH object 
(|Moriartv-Schieven et aflll992t iGomez et alj|1997l) : they 
suggest it could actually be a substellar Class I object. 
However, a recent spectral type d etermination of M4- 
M5 for this object (jLuhmanl 120061 ) suggests that it will 
probably evolve into a very low-mass star. 

Since this Class I object displays a silicate emission fea- 
ture, and in the near-IR the central binary is detected, 
its envelope likely has low densities. Our model suggests 
a reference density of 4.0 x 10" 15 g cm -3 , which is the 
lowest pi value among the protostars modeled here, indi- 
cating an advanced stage of evolution for this protostar. 
Our model further includes a relatively small centrifu- 
gal radius of 30 AU, a cavity semi-opening angle of 15°, 
and an inclination angle of 70° (see Figure [T4")) . IRAS 
04248+2612 is similar to CoKu Tau/1, which is also seen 
close to edge-on through a low-density envelope. The 
short-wavelength emission is fitted very well; it is dom- 
inated by the stellar component, which agrees with the 
fact that the central binary is detected at near-IR wave- 
lengths. The emission from far-IR to mm wavelengths 
is underestimated; however, we found no set of param- 
eters that could reproduce both the mid-IR and long- 
wavelength parts of the SED. Thus, it is likely that the 
sub-mm and mm emission is generated by an additional 
component, probably by dust in a very extended region 
of co ld dust outside the infall zone (|Jayawardhana et al.l 
12001 . 

Compared to modeling results by KCH93, who ob- 
tained two models of IRAS 04248+2612, one with 
pi = 1.0 x 10~ 13 g cm -3 and i? c =300 AU, and one with 
pi=3.2 x 10" 14 g cm" 3 and R c =70 AU, our model has 
a clearly lower density (by a factor of 8) and a smaller 
centrifuga l radius than th e lowe r values of KCH93. Also 
models by W hitney et al.l (|1997f ) resulted in a larger cen- 
trifugal radius (100 AU) and reference density (1.5 x 
10" 14 g cm" 3 ) than our values; however, they derived a 
large inclination ang le, as we did with o ur model. Recent 
modeling efforts by IStark et alj (|2006l ) yielded a model 
very similar to ours. 

Based on our fit to the IRS spectrum, especially the 
presence of a silicate emission feature, we can exclude 
a much higher density, and the fit to the shape of the 
IRS spectrum between 15 and 25 /im also constrains 
the centrifugal radius to lie around 30 AU. Compara- 
ble model fits could be achieved with p\ values in the 
(3.0-5.0) xlO" 15 g cm" 3 range and R c values between 
20 and 40 AU, with smaller reference densities accompa- 
nied by smaller centrifugal radii. Our inclination angle of 
70° not onl y matches two previously published mod els of 
this object (jWhitnev et al.lll997t IStark et al.ll2006ft . but 
is also in go od agreement with t he near-IR morphology 
observed bv lPadgett et all l| 1999ft . 

4.2.10. 04264 +2433 

The IRS spectrum of this IRAS source shows a promi- 
nent silicate emission feature, but also possibly some ice 
features in the 5-8 /jm wavelength region and an overall 



rising SED in the mid-IR. Therefore, it is probably a pro- 
tostar surrounded by envelope material. This notion is 
further supported by the fact that IRAS 04264+2433 is 
surrounded by a reflecti on nebulosity and l ikely emitting 
a small bipolar HH jet (|Devine 

Similar to IRAS 04248+2612 and CoKu Tau/1, the 
detection of a silicate emission feature in this protostar 
indicates that the envelope has low densities and is thus 
likely at an advanced evolutionary stage. We can repro- 
duce the SED of IRAS 04264+2433 fairly well with a 
model whose reference density is 5.0 x 10 -15 g cm -3 ; in 
addition, the shape of the IRS spectrum between about 
13 and 20 (im constrains the value of the centrifugal ra- 
dius, found to be 30 AU (see Figure [T5)) . The inclination 
angle of 87° and cavity semi-opening angle of 13° were 
fine-tuned to yield a good fit over the 5-20 pm. range. The 
model overestimates the emission in the near-IR some- 
what, which can be mostly be attributed to the envelope 
component. 

A previous model of IRAS 04264+2433 by KCH93 indi- 
cated a density twice as high as our value, a centrifugal 
radius of only 10 AU, and an inclination angle of 60°. 
Despite the differences, the earlier model already gauged 
the properties of this object, namely a low-density, highly 
inclined e nvelope with a small d isk. On the other hand, a 
model by Wh itney et al.l (|1997ft yielded different model 
parameters: a centrifugal radius of 50 AU, a reference 
density of 3.75 x 10~ 14 g cm" 3 , and an inclination angle 
of 32°; only their cavity semi-opening angle is similar to 
our value. Our model would yield a comparable fit with 
a centrifugal radius of 20-40 AU, but not with smaller or 
larger values for p\. 

In order to get enough flux from the far-IR to the mm 
wavelength range, we had to increase the luminosity by a 
factor of 2 compared to the bolometric luminosity of 0.37 
L Q measured by KH95. On the other hand, integrating 
the measured fluxes of the SED yields a bolometric lumi- 
nosity of 0.5 L ; given the high inclination angle of this 
object, the luminosity determined in this way is likely an 
underestimate of the true value. 

4.2.11. 04295+2251 

This object, also known as LI 536 IRS, is relatively 
faint; based on its bolometric temperature of 270 K, it is 
classified as a Class I object, but its large s ub-mm lumi- 
nosit y would suggest it is a Class object ()Young et al.l 
I2003D . Given that IRAS 04295 +2251 seemed no t to be 
extended in the near-IR or mm, iPark fc Kenvonl {2002) 
indicated that it could be an edge -on disk; however, 
data obtained by lEisner et al.l (|2005f ) shows that IRAS 
04295+2251 is somewhat extended at 1.3 mm and also 
at 0.9 jum, where scatter e d ligh t is detected. In addition, 
iMoriartv-Schieven et al.l (|1992f ) detected an outflow from 
this source, implying an early evolutionary stage for this 
object. 

Our envelope model for IRAS 04295+2251 reproduces 
the general shape of the SED (see Figure [16]), but it 
does not represent a good fit to the 10 silicate fea- 
ture and somewhat underestimates the near-IR and mm 
emission. However, at the longer wavelengths, our model 
reproduces those fluxes measured with smaller apertures, 
10"-20"in size (see Table [2]); this suggests that a region 
of cold dust outside of the infall region could be present. 

A luminosity of 0.8 Lq was necessary to obtain more 



12 



Furlan et al. 



emission especially at longer wavelengths, even though 
the source's bolometric luminosity is 0.5 L© (determined 
from integrating the fluxes of the SED). In the litera- 
ture, Lboi values v ary from 0.44 L Q (KH95) to 0.64 L 
(jMvers et al.lll987f ). A small centrifugal radius of 20 AU 
was adopted to reproduce the steep rise of the SED past 
15 /im; a value larger than 30 AU would result in a flatter 
SED in this wavelength region, while a smaller R c would 
decrease the near-IR emission and cause an even larger 
overestimate of the far-IR flux. A reference density up 
to 30% larger than our adopted value of 8.0 x 10" 15 g 
cm -3 would still yield an acceptable fit, but require lower 
inclination angles (40°-50°) to generate only a shallow 
silicate absorption at 10 /im; these models would over- 
estimate the emission between 3 and 10 pro. and in the 
far-IR. 

A previous model of IRAS 04295+2251 by KCH93 
resulted in a larger centrifugal radius, 70 AU, and a 
smaller inclination angle, i =30°, but in a compa- 
rable reference density. We can exclude an inclina- 
tion angle smaller than about 60° with the current 
model parameters, since it would increase the emission 
in the 2-8 pm range and at far-IR wavelengths, where 
the model already ov erestimates the observed emission. 
I Whitney et al.l (|1997ft modeled this object with very dif- 
ferent model parameters, including a larger centrifugal 
radius and reference de n sity, a nd an inclination angle of 
only 18°. lEisner et all (|2005l) used the same value for 
pi , 3.75 x 10~ 14 g cm" 3 , and a similar inclination angle 
as lWhitnev et al.l (|1997ft , but an R c of 30 AU, which is 
comparable to our result. 

The relatively low reference density, combined with the 
inclination of 70° to the line of sight and a cavity semi- 
opening angle of 5°, results in a self-absorbed silicate 
feature at 10 /im, with part absorption, part emission 
characteristics. However, our model does not reproduce 
the detailed shape of the observed 10 /im feature. This 
is likely due to the fact that the disk component used 
in the model does not include the optically thin disk at- 
mosphere that would generate a 10 /im silicate emission 
feature, which would then be subject to absorption by 
dust in the envelope. In the current model, the disk 
component only contributes a silicate absorption feature 
(i.e., continuum extinguished by the envelope), while the 
dust in the envelope is responsible for the emission com- 
ponent at 10 /im (see the dash-dotted and long-dashed 
lines, respectively, in Figure fTS)) . Also the shape of the 
inner cavity could affect the depth of the silicate feature, 
for example by resulting in a hotter inner region for a 
given p\. It is likely that only a more complex model 
will yield a better fit, considering that this object has a 
relatively low-density envelope and compact morphology, 
but large sub-mm luminosity. 

4.2.12. 04302+2247 

In the near-IR, this object is surrounded by two 
bright, symmetric n ebulae separated by a dark lane 
(|Padgett et alj 11999ft . Models of the very symmetric 
quadrupolar structure suggest that IRAS 04302+2247 
is s een almost exactly at an inclination angle of 
90° (jLucas fc Rochdfl997l IWolf. Padgett. & Stapelfeldtl 
2003). The central source is likely obscured by a circum- 
stellar disk, and the near-IR nebulosity, seen in scattered 
light, traces the walls of a bipolar cavity in the enve- 



lope ((Lucas fc Rochdll997l: IWolf. Padgett, fe Stapelfeldtl 
2003). Outflows oriented perpe ndicular to the dark lane 
were detected from this obj ect (jMoriartv-Schieven et al.l 
Il992t iLucas fc Roche]|1998f) . The edge-on orientation of 
IRAS 04302+2247 may explain why it is invisible in the 
optical and very faint in the mid-IR. The SED of this 
object peaks close to 100 /im and shows a pronounced 10 
/im absorption feature. 

The SED of IRAS 04302+2247 can be fitted with an 
envelope model with an inclination angle of 89°, moder- 
ate density [pi = 3.0 x 10" 14 g cm" 3 ), and a centrifugal 
radius of 300 AU (Figure [17)) . The deep silicate absorp- 
tion feature and the steep slope of the IRS spectrum be- 
yond 10 /im require a high inclination angle, consistent 
with the near-IR images. A higher density would also 
increase the depth of the silicate feature, but it would 
decrease flux levels over the 2-50 pm region, too. The 
semi-opening angle of the cavity is relativ ely wide (22°) 
also ro ughly consistent with the results of iPadgett et al.l 
(1999). The shape of the IRS spectrum beyond 12 /im, 
and the width of the SED, constrain the value of R c , 
which has to be large. Previous models by KCH 93, 
which had much fewer constraints, derived a centrifugal 
radius of only 70 AU, but a density about a factor of 3 
higher than in our mo del, and an inclination angle of 60°. 
iWhitnev et al.l (|1997f ) found an even smaller value for R c , 
10 AU, and a reference density half as large as our value, 
but their results for cavity and inclinati o n ang les agree 
with ours. Recent models bv lStark et alj (|2006ft resulted 
in very similar parameters as for our model, except for a 
somewhat smaller value for p\. 

The model emission is somewhat above the IRAS data 
points at 60 and 100 /im; this is likely due to the fact that 
our adopted luminosity of 1 L is higher than the ob- 
served bolometric luminosity of 0.34 L Q (KH95), which 
was determined from fewer flux measurements (A < 100 
/im). Our model indicates that our adopted luminosity 
might be slightly higher than the actual source luminos- 
ity. 

4.2.13. 04325+2402 

This protostar is surrounded by a bipolar reflection 
nebula that is seen in scattered light in the near-IR an d 
shows a complex morphology (Har tmann et al.l fl999). 
It is likely a multiple system, consisting of a central, 
subarcsecond binary and a faint companion about 8" 
away; both the binary and the companion are surrounded 
by an a ccretion disk and are accreting from their en- 
velopes (jHartmann et alj Il999t ) . The subarcsecond bi- 
nary, which instead could just be the result of scattering 
off a complex structure around a single star, is the likely 
source of the reflection nebula detected in the near-IR 
(jHartmann et al.l Il999f ) . This object is also associated 
with a molecular outflow (|Moriartv-Schieven et al.lll992l ; 
iHogerheiide et al.lll998l ). 

Our IRS observations were centered on the central 
source (see Figure [18|) . Even though the SL slit was 
oriented such that it also included the companion about 
8" away, this component did not enter our extraction 
window, and we also did not detect it separately; this 
is expected, since it i s faint and likely seen edge-on 
(jHartmann et al.l [19991 . and even the spectrum of the 
bright central source is only about 40 mJy between 5 
and 10 /im, with a steeply rising continuum beyond. In 
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LL, which is oriented close to perpendicular to SL, only 
the central source was included. 

An envelope model that fits this object well requires 
a sheet-collapse model (77 = 1.0), moderate density 
(pi = 3.0 x 10~ 14 g cm -3 ), a centrifugal radius of 100 
AU, a cavity semi-opening angle of 15°, and an inclina- 
tion angle of 80° (Figure [TH]) . Choosing a sheet-collapse 
model with J7 = 1.0 was necessary to fit the IRS spec- 
trum in the 15-30 fim range. The inclination angle is 
larger than the value of ~ 60° inferred from th e orienta- 
tion of the outflow bv lHogerheiide et alj (|1998f ). but the 
shape of the IRS spectrum requires a high-inclination 
model. A previous model by KCH93 also determined 
i = 60°, but a centrifugal radius and a reference density 
th ree times as lar g e as t he values derived here. Models 
bv lWhitnev et alJ (|1997l ) resulted in an R c value of only 
50 AU, but in a reference density only somewhat larger 
than our value, and very similar values for cavity and in- 
clination angles (11° and 81°, respectively). Compared 
to previous models, the IRS spectrum provides better 
constraints; a larger R c would result in a flatter mid- 
IR SED, while a larger density would generate a deeper 
silicate absorption feature. 

The silicate feature is not well fitted with this model; it 
seems to be in absorption on the short-wavelength side, 
but with an emission component on the long-wavelength 
side. The binary nature of IRAS 04325+2402 could 
account for the peculiar shape of the silicate feature, 
with one component possibly responsible for some sili- 
cate emission. 

4.2.14. 04361 +2547 

This Class I object, also known as TMR 1, is deeply 
embedded in its dense environment and su rrounded by a 
large reflection nebula, about 30" in size (|Tamura et al.l 
fl99ll). It is a close binary with a separation of 0'.'31 
(jTerebey et al.l 119981 ): the "companion" detected at a 
projected separation of 10" from this protostar, at the 
end of a nebulosity extending from the source to the 
south east, is likely a background star (jTerebey et al.l 
12000ft . IRAS 04361+2547 is e xtended in the near- 
IR, mid-IR and mm continuum (iMotte fe Andrei [20011 : 
iPark fc Kenvonll2002l ; iLuhman et al.l l2006h . An outflow 
has been detected from this object, but its morphology 
is not clear ly bipolar, suggesting a n inclination angle of 
about 60° (|Hogerheiide et al.lll998l ). 

The steeply rising SED of IRAS 04361+2547 indicates 
that this object is seen at a high inclination angle (Fig- 
ure [5D]). A good envelope model fit requires a sheet- 
collapse model (77 = 1.5), an inclination angle of 80°, 
a cavity semi-opening angle of 15°, an R c value of 100 
AU, and a fairly typical density (pi = 2.0 x 10~ 14 g 
cm -3 ). Modeling by KCH93 resulted in a reference den- 
sity of 3.2 xl0~ 14 g cm~ 3 and a centrifugal radius of 
only 10 AU. Also, while our model suggests a close to 
edge-on orientation, KCH93 derived an inclination angle 
of only 30° . La t er mo dels bv lWhitnev et all (|1997l ) and 
iGramajo et al.l (|2007f ) confirm our result of a large in- 
clination angle (and also our cavity opening angle), but 
they derive centrifugal radii about half as large as our 
value, and generally a somewhat larger reference density. 

Our adopted reference density is required to fit the 
SED beyond 50 pm, as well as to produce a silicate fea- 
ture at 10 pm that is mostly dominated by an emission 



component generated in the envelope (long-dashed line 
in Figure l20|). The R c value of 100 AU was chosen to 
reproduce the SED shape from about 15 to 50 pm. A 
larger R c would result in a flatter SED, while a smaller 
value would result in more mid-IR flux and a pronounced 
silicate absorption feature. Choosing Tj s t ar = 0.9 and 9 = 
15°, combined with the high inclination angle, depresses 
the emission below 10 pm. Our model does not repro- 
duce the deep ice features at 6.0 and 6.8 pm, but it traces 
the continuum level between 5 and 8 pm (where emission 
from the disk dominates). We note that the IRS fluxes 
are lower than the IRAC data points; this is likely due to 
the e xtended emission detected in IRAC (jLuhman et al.l 
|2006|) and thus an aperture effect. Below 2 pm, the near- 
IR scattered light component from the envelope is likely 
overestimated. We note that even though the luminosity 
of our model, 4.0 L Q , is larger than the measured bolo- 
metric luminosity of 2.5 L Q (see Table [4]), it is necessary 
to yield enough flux at all wavelengths. 

4.2.15. 04365+2535 

This IRAS source is also known as TMC 1A; it is op- 
tically invisible and thus deeply embedded. It is sur- 
rounded by a cometary ne bulosity in the near-infrared 
(|Tamura et al.l 119911 11996D , and it is e xtended both in 
the near-IR and at mm wavelengt hs (IPark fc Kenvon 
2002). Both IChandler et all (|1996h and iTamura et a! 
(1996) measured a bipolar outflow which is somewhat 
conical in shape due to the blueshifted lobe being more 
prominent than the redshifted one. The near-IR nebulos- 
ity is likely light from the central sou rce scattered by dust 
in the outflows (jTamura et al.|[l9 96). Also, dust contin- 
uum measured at 790 /im and 1.1 mm is aligned with the 
outflo w direction ([Chandler et al.l[l~998h . IChandler et al.l 
( 1996) concluded that the inclination angle of this object 
lies between 40° and 68°, with a cavity of semi-opening 
angle between 15° and 21°. 

The SED of 04365+2535 is fitted by an envelope model 
with relatively high density (pi = 4.5 x 10~ 14 g cm" 3 ), 
a centrifugal radius of 50 AU, a small cavity, and an in- 
clination angle of 30° (Figure l2Tj) . Compared to models 
by KCH93, both our reference density and centrifugal 
radius are smaller, by factors of 0.45 and 0.17, respec- 
tively, and our model suggests a more face-on orientation 
as opposed to the more highly inclined model (i = 60°) 
of KCH93 and the inclination a ngle determined from the 
outflows (IChandler et al.lll996f). More recent m odels by 
IWhitnev et al.l([1997l ) and lGramaio et alJ (|2007h resulted 
in a large inclination angle (70°-80°), but in values for 
R c and p\ very comparable to ours. 

The high density is required to fit the deep silicate ab- 
sorption feature, as well as to shift the peak of the SED 
to longer wavelengths. Our model does not fit the sub- 
mm and mm fluxes measured with a large aperture (40"- 
60"), but matches the measurements done with apertures 
< 20" (see Tabled]). This suggest the presence of some 
extended, cold dust beyond the infalling envelope. A cen- 
trifugal radius of 50 AU was chosen to fit the emission 
from the silicate absorption feature out to about 50 /im; 
a model with larger values for p\ and R c (~ 6.0 x 10 14 
g cm~ 3 and 70 AU, respectively) would also yield a good 
fit, but underestimate the emission in the 20-30 pm range 
somewhat. An even larger centrifugal radius would de- 
crease the emission from 15 to 40 pm. The inclination 
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angle of o ur envelope model is sm aller than the value sug- 
gested bv lChandler et all (|1996h and by previous models; 
however, we can exclude a larger inclination angle, since 
it would result in a steeper SED from 15 to 60 yum, as well 
as less flux in the 2 to 10 /mi region. We can also exclude 
a cavity semi-opening angle larger than about 5° , since it 
would decrease the flux in the 5-15 /im wavelength range. 

4.2.16. 04368+2557 

This object has a bolometric temperature of less than 
70 K, thus placing it in the Class d omain, the earlier , 
deeply embedded pro tostellar stage (|Chen et al.l 119951 : 
iMotte fc And re 2001). It is very fa int and extended in 
the near-IR (|Park fe Kenvonll2002l ): no point source is 
detected in 2 MASS. The if- b and d ata point in our SED 
plot is from IWhitnev et alJ (|1997l ). who only detected 
some faint reflection nebulosities at the position of the 
IRAS source. Since the cent r al sou rce is not detected 
in the near-IR, iTamura et al.l (p~996) conclude that this 
object might be seen edge-on, obscured by a disk. 

Our IRS observations were centered on the position 
de termined by i n terfero metric observations at 2.7 mm 
by lOhashi et al.l (|1997af ) , which also coincides with the 
po sition of the source in high-resolution 7 mm VLA maps 
bv lLoinard et al.l {2002). The latter authors resolved this 
IRAS source into two components separated by 0'.'17 and 
suggested that it is a binary in which one component 
is surrounded by an almost edge-on disk. In addition, 
this object might have a faint companion located 20" 
to the northwest of the main source, but it was only 
detected at 800 /mi ([Fuller. Ladd. fc HodaDp|[T9 96) and 
not c onfirmed at 1.3 mmjMotte fc Andrell2001l) "and 2.7 
mm (jOhashi et al.lll997af) . 

IRAS 04368+2557 drives molecular outflows, which 
are located symmetrically around the central source and 
also perpendicular t o the elongated source det ected by 
lLoinard et all (|2002l ). iHogerheijde et all (|1998| ) inferred 
an inc lination angle i > 65° from observations of the out- 
flows. lOhashi et al.l (|1997af ) observed this object at high 
resolution in the millimeter wavelength range and con- 
cluded that it is surrounded by an elongated envelope 
~ 2000 AU in radius with large bipolar cavities, seen 
edge-on. 

The SED of IRAS 04368+2557 (Figure[22]) is similar to 
that of IRAS 04302+2247 (FigureQH); both objects have 
a deep silicate absorption feature and SEDs that peak 
between 60 and 70 /im. Both are seen close to edge- 
on, but the reference density of the IRAS 04368+2557 
model fit is higher (pi = 4.0 x 10~ 14 g cm~ 3 ), and the 
centrifugal radius smaller (R c = 200 AU). On the other 
hand, about 20% changes in p\ and R c yield comparable 
model fits, as long as the cavity semi-opening angle lies 
between 20° and 30°, and the source is seen edge-on. We 
can exclude a high density as derived by KCH93 (3.2 
xlO -13 g cm -3 ), but agree reasonably well with their 
result of R c (300 AU) and inclination angle (60°-90°). 

In the near-IR, our model reproduces the K-band mea- 
surement, but seems to overesti mate the 3.6 and 4.5 /im 
emission, as measured by IRAC (Hartma nn et al.l l2005). 
However, this object is extended in IRAC images, and 
the photometry from the 6" aperture shown in Figure [22] 
doe s not capture all of the strong extended emission (see 
also lHartmann et al.l 120051 ) . 

As with IRAS 04302+2247, the values for p x and R c 



are fairly well constrained by the depth of the silicate ab- 
sorption feature, the shape of the IRS spectrum beyond 
10 /im, and the width and peak of the SED (the 100, 
160, and 1300 /im measurements are likely overestimates 
due to aperture effects). The similarity between the two 
IRAS sources suggests that some Class I objects, seen at 
high inclination angles, could resemble Class objects. 

4.2.17. 04381+2540 

This IRAS source is also known as TMC 1; it 
is in visible in the optical and thu s deeply embedded 
(e.g.. iMoriartv-Schieven et al.l ll~992). It is surrounded 
by a near-infrared n ebulosity seen in s cattered light 
(jTamura et al.l 119911: ITerebev et all 120061 ). The coni- 
cal shape of this near-infrared emission indicates that 
it likely traces an evacuated outflow cavity; the CO 
molecular outflow has a semi-opening angle between 13° 
and 19° and an inclin ation angle between 40° and 70° 
([Chandler et al.lfl99ll . 

Dust emission orie nted perpendicular to outflows was 
detected at 790 /im ([Chandler et al.lll998[ ). lYoung et all 
(|2003| ) measured extended 450 and 850 /xm emission from 
this object; based on the source's properties and model 
calculations, they confirmed that this object is a Class I 
sourc e. From radial i n tensit y profile fits to the 1.3 mm 
map, IMotte fc Andrei (|200fl) derived an outer envelope 
radius of 3900 AU. Recently, lApai et all (|2005f) found a 
0'.'6 companion to IRAS 04381+2540; they suggest that 
this companion might be a young brown dwarf. 

A good model fit for IRAS 04381+2540 (see Figure 
I23p requires a somewhat higher luminosity than its bolo- 
metric luminosity of 0.7 L (see Table [4j to be able to 
reproduce the long-wavelength emission. By increasing 
the inner disk radius to 5 stellar radii, the emission in the 
near-IR and mid-IR out to 8 /im is decreased, resulting in 
a better fit. The centrifugal radius of 70 AU and the ref- 
erence density of 3.0 x 10~ 14 g cm~ 3 are well constrained 
by the depth of the silicate absorption feature and the 
shape of the IRS spectrum from 15 to 30 /im. Decreas- 
ing R c to 60 AU and the inclination angle to 35° yields 
a nearly identical fit, but we can exclude larger changes 
in these parameters, as well for the reference density. 

Previous modeling by KCH93 yielded a value of 1.0 x 
10~ 13 g cm~ 3 for pi and 300 AU for R c , and an incli- 
na tion angle of 30° . How ever, better-constrained models 
by IWhitnev et al.l (I1997D resulted in values for p\ and 
R c similar to ours. lEisner et al.l (|2005f ) derived a smaller 
centrifugal radius (30 AU) and also a larger reference 
density ( 6.75 x 10~ 14 g cm" 3 ). 

Our results for cavity semi-opening angle (0 =10°) and 
inclination (i =40°) are at th e lower limits of the range 
suggested by IChandler et al.l ((1996). O ur value for 9 
agrees well with the modeling result of IWhitnev et al.l 
(fl997h . but our inclination angle is smaller; on the other 
h and, this latter quant ity is consistent with the result 
of ITerebev et al.l ([20061 ). but not with their cavity semi- 
opening angle of 40°. Higher values for both parameters 
can be excluded based on the constraints placed by our 
IRS spectrum on the envelope models. Finally, we note 
that we adopted an outer envelope radi us of 10000 AU, 
which is larger than the value derived bv IMotte fc Andra 
(2001), but necessary to fit the long- wavelength part of 
the SED. 
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4.2.18. 04489+3042 

Based on the M6 spectral type of IRAS 04489+3042, 
iWhite fc Hillenbrand] (|2004D indicated that it could be 
substellar, a Class I brown dw arf. How e ver, recent 
spectral type determinations by iLuhmanl ([2006) sug- 
gest that its spectral type is M3-M4, thus placing it 
into the low- mass star regime. IR AS 04489+3042 is not 
extended in the near-IR or mm (|Park fc Kenvonl [2002; 
iMotte fe Andre1l200lD, and no outflow has been detected 
(jGomez et al.lll997| ). iPark fc Kenvonl (|2002f ) suggest it 
could be in transition between Class I and II stage, prob- 
ably seen close to pole-on. 

Our preliminary envelope model fit for this object (Fig- 
urc[24|) suggests a low density (pi = 1.0 x 10~ 14 g cm -3 ), 
a small centrifugal radius (R c — 15 AU), a very small cav- 
ity, and a low inclination. Models by KCH93 yielded the 
same reference density as our model and a similar low in- 
clination angle, but a centrifuga l radius of 70 AU. On the 
other hand, modeling efforts bv lWhitnev et all (| 19971 ) re- 
sulted in larger values for all these three parameters. We 
can exclude a centrifugal radius larger than about 25 AU, 
since it would decrease the flux over the 10-60 /jm spec- 
tral range. The low density and inclination angle of our 
model are mainly set by the very weak silicate absorp- 
tion feature and roughly flat SED from about 5 to 15 
/im. The model does not yield a good fit of the details 
of the IRS spectrum below about 15 /zm, but the overall 
shape of the spectrum is reproduced reasonably well. To 
fit the data point at 1.3 mm with the set of model pa- 
rameters described above, we adopted an outer envelope 
radius of only 1000 AU; the small envelope contributes 
only little emission in the mm, where the emission from 
the disk dominates (see the dash-dotted line in Figure 
I24j) . Since there are no sub- mm measurements available, 
and the IRS spectrum is relatively flat with no strong 
silicate feature, the current model parameters are only 
weakly constrained. 

4.2.19. CoKu Tau/1 

Coku Tau/1 was identified as a Class II object by 
KH95, but it is surrounded by filamentary reflection neb- 
ulae detected in the near-IR and suggested to rep resent 
the walls of outflow cavities ([Padgett et al.l fl999). It is 
seen close to edge-on, and, like IRAS 042 48+2612, is 
a subarcsecond binary dPadgett et al. I [19991 ). Similar to 
IRAS 04248+2612 and 04264+2433, CoKu Tau/1 has a 
silicate emission feature and an infrared SED that is ris- 
ing from 3 to about 20 /*m (see Figure [3]). 

The appearance of CoKu Tau/1 in the near-IR suggests 
that it is surrounded by a low-density envelope and an 
edge-on disk. It is likely a more evolved protostar, at 
a stage when the envelope has begun dissipating. The 
binary and the inner disk arc extinguished by the outer 
disk, but since in the near-IR the optical depth of the 
outer regions of the circumbinary disk decreases, the two 
stars at the center are seen at near-IR wavelengths. As in 
the optical, there is likely also some stellar light scattered 
by the envelope contributing to the near-IR emission. 

Our model of CoKu Tau/1 (see Figure |2"5]) requires a 
low density (pi = 5.0 x 10~ 15 g cm -3 ), a small centrifu- 
gal radius (R c — 40 AU), a cavity semi-opening angle of 
5°, and a high inclination angle (i = 80°). We adopted 
a luminosity of 1.1 L Q , which is slightly higher than the 



measured bolometric luminosity of 1.0 L Q (see Table [4]). 
The inner disk radius was set at 7 stellar radii to de- 
crease the emission from the near-IR to 8 /im, which is 
still higher than the observed fluxes, most likely due to 
an overestimate of the scattered light emission from the 
envelope. 

The shape of the IRS spectrum between 12 and 20 /im 
constrains the centrifugal radius, while the low density is 
required to generate a silicate emission feature and cause 
the SED to peak around 30 /im. Decreasing p\ by about 
20% and increasing or decreasing R c by roughly the same 
percentage yields comparable fits, provided the inclina- 
tion angle lies in t he 70°-80° range. A previous model by 
IStark et"aLl (|2006[ ) yielded a reference density smaller by 
about a factor of three, a comparable values for i? c , but 
a somewhat smaller inclination angle and a larger cavity 
opening angle. We can exclude a larger cavity, since it 
would decrease the flux of the 10 /jm silicate emission 
feature. 

As mentioned earlier, the silicate emission feature at 
10 ura. is generated by the envelope. The disk just con- 
tributes a weak silicate absorption feature, which results 
from our simple disk model that includes only continuum 
emission, extinguished by the envelope. A more realis- 
tic edge-on disk would yield a somewhat stronger ab- 
sorption feature; however, the envelope emission clearly 
dominates over that of the disk at 10 /im, a fact that 
also applies to IRAS 04264+2433 and, to a lesser extent, 
IRAS 04248+2612 (long-dashed gray lines in Figures QH 
[l"5" l [25 ]) . Thus, the silicate emission feature is not greatly 
affected by the shape of the silicate feature originating 
from the disk, and our current models of low-density en- 
velopes have some validity. 

4.2.20. DG Tau B 

This object displays some similarity to HH 30; its near- 
IR HST images show an edge-on disk and a bipolar re- 
flection nebulosity (see Fig. [I]), which was interprete d as 
the walls of an outflow cavity ([Padgett et al.l [l999h . In 
the optic al, the nebula is visible, too, but the star is not 
detected ([Stapelfeldt et al.lll997n . DG Tau B is also the 
source of a jet, which is oriented perpendicular to the disk 
and along the axis of sy mmetry of the reflection nebula 
([Eisloffel fc Mundtlll998| ). A molecular outflow, oriented 
along the redshifted optical jet, has also been observed 
(|Mitchell. Sargent, fc Manningslll997ft . 

Our preliminary model for DG Tau B (see Figure I26p 
has an inclination angle of 55° and a cavity semi-opening 
angle of 10°; both values are smaller than expected from 
the optical and near-IR images. However, a more highly 
inclined envelope would result in a steeply rising SED 
over the mid-IR range, which would not match the ob- 
served spectrum. The centrifugal radius (and thus outer 
disk radius) of 60 AU is poorly constrained; a larger value 
would better fit the slope of the IRS spectrum beyond 
20 /itm, but result in a shallower silicate feature at 10 
/im. The deep silicate absorption feature requires a rel- 
atively high density (pi = 3.5 x 10~ 14 g cm -3 ), but its 
long-wavelength wing is not fitted very well. In addition, 
the model is not well-constrained at longer wavelengths, 
since DG Tau B lacks any sub-mm measurements, and 
the IRAS fluxes at 60 and 100 /im are uncertain due to 
the proximity of DG Tau: the far-IR fluxes contain the 
emission of both DG Tau B and DG Tau, which lies 1' 
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to the northeast of DG Tau B and is about two to three 
times brighter in th e mid-IR. 

Recent models bv lStark et al.l (|2006l ) yielded quite dif- 
ferent model parameters, except for the reference den- 
sity; they were able to fit the near-IR images of DG Tau 
B adopting a luminosity of 0.2 L©, a centrifugal radius 
of 300 AU, a pi value of 3.75 x 10~ 14 g cm -3 , and cav- 
ity semi-opening and inclination angles of 30° and 73°, 
respectively. This set of parameters would not fit the 
IRS spectrum. In particular, we adopted a luminosity of 
2.5 L©; a smaller luminosity, closer to the value derived 
from integrating under the SED (1.8 L©), would result in 
less far-IR flux and require a smaller inclination angle to 
still reproduce the mid-infrared spectrum, which would 
be inconsistent with the HST images, which indicate a 
more edge-on orientation. 

4.2.21. HL Tau 

Even though HL Tau is considered by some a clas- 
sical T Tauri star surrounded by an a ccretion disk, it 
is also embedded in an envelope (e.g., iMotte fc Andrei 
120011: 1 White fc Hillenbrandll200l : this is supported by 
our IRS spectrum, which shows the silicate feature in 
absorption, as well as ice absorption features in the 5-8 
pm range and at 15.2 pm. HL Tau is also very bright 
at sub-mm and mm wavelengths, and a rotating, elon- 
gated, ~ 2000 AU long structure has been mapped i n 
13 CO dBeckwith et all 11995 iSareent fc Beckwithlfl99l . 
iStanelfeldt et all (|1995lT detected a compact, one-sided 
reflection nebulosity in the optical; the star is not de- 
tected. HL Tau drives a powerful, extend ed molecular 
outflow which was mappe d in CO emission onin et al.l 
I1996D . IClose et al.l (| 19971) observed HL Tau in the near- 
IR with adaptive optics and resolved the inner accretion 
disk and bipolar outflow cavities. They concluded that 
HL Tau is surrounded by an accretion disk of about 150 
AU radius, inclined at 67°, and an infalling envelope of 
about 1200 AU radius. However, the environment of HL 
Tau is quite complicated, partly due to its proximity to 
XZ Tau, which is likely at the origin of an ex panding shell 
that compresses the nebula around HL Tau (Welc h et alJ 
l2000l) . 

The SED of HL Tau can be reproduced with an enve- 
lope model with a luminosity of 8.0 L©, which is some- 
what larger than its measured bolometric luminosity (see 
Table |U), a relatively high density (pi — 4.5 x 10 -14 g 
cm -3 ), an R c value of 100 AU, and flattening parameter 
Tj of 1.0, i.e., a sheet-collapse model (Figure 127]) . The 
latter quantity was adjusted to yield more flux in the 12 
to 30 pm region compared to the TSC models. A lower 
density would yield a 10 pm silicate feature that is too 
narrow, while a density up to 6.0 x 10~ 14 g cm -3 would 
still yield a good fit. An even higher density would re- 
quire a larger centrifugal radius, which would result in 
an SED that overestimates the flux from about 25 to 200 
pm. The R c value is constrained by the shape of the 
IRS spectrum; we can rule out R c values smaller than 
100 AU and larger than about 200 AU. These results are 
also roughly c onsist ent with previous modeling results by 
ICalvet et alJ ((1994), who adopted a reference density of 
3.2 xlO -14 g cm -3 and a centrifugal radius of 200 AU. 

We note that, despite the finding bv I Close etal] ([19971 ) 
that the inclination angle of HL Tau is 67° and that large 
cavities are present, our model indicates virtually no cav- 



ity an d a low inclination. The model by ICalvet et al.l 
(|1994| ) also suggests an outflow hole with a semi-opening 
angle of 10°; including such a cavity in our model would 
reduce the flux over the IRS wavelength range, but also 
increase the flux somewhat in the near-IR, where our 
model underestimates the emission. A larger inclina- 
tion angle would decrease the short-wavelength flux even 
more and generate a deeper silicate absorption feature. 
On the other hand, our value for the incl ination angle 
is clos er to the value of 30° determined bv ICalvet et al.l 
(1994), and our disk radius (~ R c ) is also r oughly con- 
sistent with the value of 150 AU found by IClose et al.l 
(fl99l . 

Our model underestimates the emission at sub- 
mm and mm wavelengths, despite an outer enve- 
lope radius of 10,000 AU. Models of HL Tau by 
iD'Alessio. Calvet. fc Hartmannl (|1997l ) have shown that 
the outer disk requires higher temperatures to reproduce 
the SED shape at longer wavelengths, implying irradia- 
tion of the circumstellar disk by the infalling envelope. 
Since our model does not include this additional heat- 
ing, our disk emission is likely underestimated; thus, the 
"missing" flux at the longer wavelengths could be at- 
tributed in part to this additional disk component. In 
addition, dust outside the infall region could contribute 
to the mm emission. 

4.2.22. L1551 IRS 5 

This object is the most luminous Class I object in Tau- 
rus (Lboi = 22-28 L©; see TableSJ and is therefore a well- 
studied protostellar system. It consists of a 0'.'3 binary, 
whose components are each surrounded by a circumstel- 
lar disk , and of a circumbinary disk and a flattened en- 
velop e (jLoonev. Mundv. fc Welcblfl997t iRodriguez et all 
|1998| ). It is surrounded by an extende d reflection neb- 
ula a t optical and near-IR wavelengths ([Campbell et al.l 
119881: iTamura et all Il991t IWhite et all l2000f ). and h 
powers a highly collim ated bipolar molecular o utflow 
(|Snell fc Schioerbl fl985l : iMoriartv-Schieven et al.lfl992ll . 
Even though no optical outburst has been observed, 
L1551 IRS 5 is considered an FU Ori object, implying 
episodic events of high disk accretion and outflow rates 
([Hartmann fc Kenvonlll996D . 

The mid- and far-IR spectrum ( 2-200 pm) of this 
bright Class I object was obtained bv lWhite et al.l (|2000l ) 
using ISO; they detected gas-phase emission lines and ice 
absorption features. We also detect strong ice absorption 
features in our IRS spectrum, in addition to a very deep 
and broad silicate absorption feature (see Figure [2b) . 

L1551 IR S 5 has been modeled in detail by 
lOsorio et al.l ([2003), taking into account emission from 
the circumstellar and circumbinary disks, as well as the 
envelope. The model presented belo w does not include al l 
the components of the system as in lOsorio et al.l (|2003h 
and should therefore be considered as preliminary. 

L1551 IRS 5 is the most luminous object in our sample 
(adopted luminosity of 25 L©). Its deep silicate absorp- 
tion feature requires a high density (pi — 7.0 x 10~ 14 
g cm -3 ), and the steep SED beyond 15 pm a centrifu- 
gal radius of 100 AU (see Figure [55]). An R c value 
larger than about 150 AU would result in a decrease in 
flux between 15 and 30 pm and thus degrade the fit to 
the IRS spectrum. Our model indicates that the cav- 
ity is small (9 =5°), and the object is seen at an inter- 
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mediate inclination angle. Increasing p\ would require 
a smaller inclination angle, but a larger inclination is 
more in accordance wit h the observed outflows (i=65°; 
iHogerheiide et al.lll998t ). 

Compared to the models by KCH93, which suggested a 
reference density of ~ 10~ 13 g cm -3 , a centrifugal radius 
between 70 and 300 AU, and an inclination angle between 
30° and 60°, our model has a smaller, but still high, ref- 
erence density, an R c value towards the lower end of the 
KCH93 range, and an inclination angle in the middle o f 
the range determin ed by KCH93. IWhitnev et alj (jl997l ) 
and, more recently, iGramajo et al.1 (|2007h . modeled this 
object with a reference density of only 3.75 x 10~ 14 g 
cm" 3 , R c ~ 40 AU, 6>=20°, and i > 70°. 

Our m odel parameters ar e different from the ones 
found by lOsorio et al.l (|2003f ) , who adopted an inclina- 
tion angle of 50°, p x = 4.0 x 10" 13 g cm" 3 , R c = 300 
AU, and a flattening parameter r\ = 2.5 for the enve- 
lope. We could obtain a comparable fit of the SED out 
to 30 pm by using L = 25 L©, pi — 2.5 x 10~ 13 g cm" 3 , 
R c = 100 AU, n = 2.0, 7] star =0.l, 0=0.1°, and i =45°, 
but the far-IR emission would be overestimated. On the 
other hand, the sub-mm and mm emission would be bet- 
ter matched with the sheet-collapse model. Thus, the 
SED of L1551 IRS 5 can be reproduced by a TSC model, 
but the long-wavelength emission suggests the presence 
of an additional, cold, likely elongated, dust component 
that is outside of the infall region. 

5. DISCUSSION 

5.1. Identification of Class I Objects 

The mid-infrared spectra of Class I objects have quite 
different appearances; while some differences are caused 
by individual envelope parameters like the density and 
size of the infall region, the inclination angle also plays an 
important role. A Class I object viewed pole-on through 
a cavity in the envelope can look like a Class II object 
with some additional long-wavelen gth excess emissio n 
(as some FU Ori objects appear; see lGreen et al.ll2~006l) . 
while an edge-on Class I object can have the appear- 
ance of a Class object. However, the envelope around 
a true protostar will generate excess emission at mid-IR 
to mm wavelengths, resulting in a larger infrared excess 
than expected from an accretion disk alone, and less ex- 
cess when compared to the large, cool envelopes around 
Class objects that peak at sub-mm to mm wavelengths. 
In addition, if seen within a narrow range of inclination 
angles, an edge-on Class II o bject can appear simi lar to 
an embedded protostar (see iD'Alessio et aT1ll999t ). but 
it would lack additional excess emission caused by the 
envelope. 

Observations at other wavelengths help in identify- 
ing the evolutionary state of an object: as noted by 
iPark &; Kenvonl (|2002D . true Class I sources are usually 
extended in the near-IR and at sub-mm or mm wave- 
lengths, the former due to light scattered by dust in in- 
ner envelope regions, the latter due to thermal emission 
by dust in the outer parts of the envelope. Furthermore, 
Class I objects usually have molecular outflows or jets 
detected in the optical/near-IR and especially at mm 
wavelengths. Objects in transition between the Class 
I and Class II stage are expected to be surrounded by 
less envelope material, in accordance with larger cavities 
in their envelopes. 



Of the 28 Class I objects presented in this paper, 7 were 
previously identified as Class II objects in KH95: IRAS 
04154+2823, 04158+2805, 04278+2253, CoKu Tau/1, 
HL Tau, IC 2087 IR, and LkHa 358. While the SED 
shape of these objects might indicate that they are T 
Tauri stars whose envelope has already mostly dissipated, 
the presence of ice features in their mid-infrared spec- 
trum, as well as the detection of extended reflection neb- 
ulosities and outflows around some of these objects, sug- 
gests that there is still some remnant envelope material. 

IRAS 04154+2823, 04278+2253, 04489+3042, IC 2087 
IR, and Lk Ha 358, as wel l as som e of the Class II objects 
presented in lFurlan et all (|2006h (IRAS 04187+1927, DG 
Tau, DP Tau, FS Tau, HN Tau, T Tau, and XZ Tau), 
could be objects in transition between the Class I and 
II stage, when a protostar has almost completely cleared 
its envelope, but is still surrounded by some nebulosity. 
In some cases, the "evolved" protostar is still generating 
outflows, but due to envelope clearing it appears more 
and more like a classical T Tauri star. In some other 
cases, like T Tau, HN Tau, and XZ Tau, a binary consists 
of a T Tauri star and a more deeply embedded source, 
which is likely an effect of the particular viewing geome- 
try of the system: either one component is more aligned 
with the outflow cavity and thus less embedded, or the 
more extinguished component is actually a T Tauri star 
oriented edge-on. Therefore, as stated earlier, the incli- 
nation angle is an important factor in determining the 
appearance of Class I and II objects. 

5.2. Ices in the Envelopes of Class I Objects 

The ice and silicate absorption features detected in the 
IRS spectra are also a useful tool to characterize Class I 
objects. The peak optical depth of these features reveals 
their strength and can therefore be used as an indication 
for their origin. To determine the optical depth of the 
CO2 ice feature at 15.2 pm and that of the silicate ab- 
sorption feature around 10 pm, we first used a spline fit 
to determine the underlying continuum. Then we derived 
the optical depth, r, assuming F b s = F c e~ T , where F b s 
and F c are the observed and the underlying continuum 
flux, respectively. The measured peak optical depths 
of the CO2 and silicate absorption features are shown 
in Figure [29], tog ether with previous m easurements by 
iGibb et alTij2004h and lAlexander et al.l (|2003l ). who de- 



termined the optical depths of silicate and ice absorp- 
tion features in ISO spectra of embedded, massive YSOs 
and low- and interm e diate- mass YSOs, respectively. For 
the lAlexander et all (120031) d ata, w e used the same con- 
version as in IWatson et all (|2004h to convert the CO2 
equivalent widths to peak opt ical depths. 

As already pointed out by IWatson et al.l (|2004h , the 
low-mass Class I objects in Taurus generally have large 
peak CO2 optical depths, especially considering that 
their peak silicate optical depth lies below ~ 2 (the ob- 
ject with the largest silicate and CO2 optical depths is 
IRAS 04368+255 7, a Class object). M ost of the ob- 
jects observed bv lAlexander et all (|2003l ) have similarly 
strong CO2 ice features for larger silicate optical depths, 
suggesting additional extinction along the line of sight 
towards these objects. It is expected that the two opti- 
cal depths track each other along lines of sight through 
molecular clouds (with some dispersion due to differ- 
ent conditions in different clouds), since ice mantles are 
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thought to grow in a similar manner onto dust grains in 
all comparably shielded regions of molecular clouds. 

The gray dashed line in Figure [55] is drawn by eye to 
roughly separate the region occupied by the Taurus Class 
I objects, where extinction from the envelope dominates 
the optical depth, from that occupied by objects seen 
along molecular cloud lines of sight, where the extinction 
is mostly due to intervening molecular cloud material. 
The boundary is not well-defined, since some of the ISO 
observations did include objects in which local extinction 
by an envelope dominates, but it separates all but one of 
the Taurus Class I objects from the majority of the ISO 
targets. 

Thu s, we confirm the conclusions of IWatson et al.l 
(2004) that dust grains in protostellar envelopes have 
larger ice mantles than grains in the ambient molecular 
cloud, implying ice mantle growth inside the envelopes. 
We defer a more detailed study of the ice absorption 
features, including the detailed composition and thus de- 
rived temperature of the ices, to a future paper (Zasowski 
et al. 2008, in preparation). 

We note that the objects with a peak silicate optical 
depth of are protostars in which the silicate feature is 
either in emission or displays both an emission and an 
absorption component. As discussed in § 13.11 the latter 
type of objects is likely seen along lines of sight that suf- 
fer from less extinction by the envelope, caused by a more 
face-on orientation or a lower-density envelope. There- 
fore, for these objects the optical depth of the silicate 
feature, taken as a measure of dust absorption in the en- 
velope, is probably underestimated due to contamination 
by the silicate emission arising from the disk. 

The one Class I object that lies in the molecular cloud 
region of the plot is DG Tau/B; it has a deep silicate, 
but only a weak CO2 ice absorption feature. This could 
indicate that grain mantles in its envelope are less rich 
in CO2 ice, but on the other hand this feature is su- 
perposed on the broad and deep silicate absorption at 
18 /im, which adds additional uncertainty to the optical 
depth measurement of the CO2 feature. 

The three objects identified as "Class I /II objects" in 
Figure[29]are IRAS 04278+2253, LkHa 358, and IC 2087 
IR. The former two objects have the lowest CO2 feature 
strengths of our Taurus sample and virtually no silicate 
absorption feature, which confirms that they are more 
evolved Class I objects. IC 2087 IR has both C0 2 and 
silicate peak optical depths comparable to those of typ- 
ical Class I objects; thus, it is still likely embedded to 
some extent. This shows that the optical depths of the 
silicate and ice absorption features can help to recognize 
more evolved Class I objects, but the effects of the en- 
vironment and viewing angle also have to be taken into 
account. 

5.3. Summary of the Models 

Envelope models aid in the identification and charac- 
terization of Class I objects; when models fit the SED 
from the near-IR to the mm, they reveal both the large- 
and small-scale structure of a protostar, its accretion disk 
and its envelope. The IRS spectrum, which covers the 
SED from 5 to 40 /im, is a particularly useful constraint, 
since the depth of the silicate feature at 10 /im is con- 
nected to the density and inclination angle of the en- 
velope, while the slope between 15 and 40 /im helps to 



constrain the envelope shape, centrifugal radius and in 
part the size of the cavity. Not all parameters can be 
uniquely defined by a model fit to the IRS spectrum, but 
together with shorter- and longer-wavelength data, the 
parameter space can be narrowed down considerably. 

We point out that, given the large number of parame- 
ters for each model, we did not determine the goodness 
of a model fit by ^-minimization, but rather by judg- 
ing by eye, based on our experience on how the different 
parameters affect the resulting model SED. However, we 
justify the choice of specific parameters and explain how 
the parameter space was narrowed. In addition, the more 
constraints a model has, in particular the availability of 
sub-mm and mm flux measurements and the source's lu- 
minosity, the more limited the choice of parameters be- 
comes. The fitting by eye is much more efficient than 
^-minimization, and in most cases we derive good fits 
that are likely minima in x 2 space. Typical uncertainties 
for our model parameters are listed in the footnote of 
Tabled 

We also note that several assumptions made in the 
modeling code could be refined to better reproduce the 
observations: e.g., the assumption that the cavity follows 
the shape of the streamlines of infalling particles, the 
assumption of a spherically symmetric density distribu- 
tion when computing the equilibrium temperature in the 
envelope, and the adoption of only the inner, optically 
thick regions of a flat accretion disk as the disk compo- 
nent. Some of our 22 models of Class I objects should be 
treated as preliminary, given the assumptions made in 
the modeling code and/or some weaker constraints due 
to the lack of data, but our models already reveal a large 
range in parameters necessary to obtain fits (see Table [3] 
for a summary of the model parameters). 

We computed models for the majority of Class I objects 
in our sample; however, we did not attempt to model ob- 
jects with poor observational constraints for their SEDs, 
or objects which are likely more evolved Class I objects 
and for which these envelope models might not apply. 
Thus, we did not attempt to model HH 30, whose spec- 
trum is very noisy and whose overall SED is poorly con- 
strained by the scarcity of measurements at other wave- 
lengths. Also IRAS 04181+2654 B, a member of a binary 
system, was not modeled, since it lacks long-wavelength 
data (we modeled component A, which was observed at 
1.3 mm). The three sources with decreasing SEDs over 
the IRS spectral range, IRAS 04278+2253, IC 2087 IR, 
and LkHa 358, were also not modeled. Finally, no model 
was produced for GV Tau, since its IRS spectrum is un- 
certain, and the fact that the two components are sepa- 
rated by l'/3 (which corresponds to about 180 AU at the 
distance of the Taurus star-forming region) likely intro- 
duces asymmetries in the envelope of the more embedded 
component. 

Our envelope models indicate centrifugal radii, and 
thus disk sizes, that cover a range from 10 AU up to 300 
AU, with a median value of 60 AU. Since the centrifugal 
radius is thought to increase rapidly with time (oc i 3 ), 
we could conclude that the objects in our sample span 
somewhat different ages (about a factor of 3). However, 
the main reason for the differences in centrifugal radii is 
likely different initial conditions of the collapse, since the 
initial centrifugal radius strongly depends on the mass of 
the central condensation and the initial angular velocity 
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of the cloud core (R c oc M 3 fl 2 ; iTerebev. Shu, fc Cassenl 
1984). In addition, some of the more deeply embedded 
protostars that have been identified as Class 0/1 objects, 
like IRAS 04166+2706 and 04368+2557, have large R c 
values, which would not be expected if the spread in 
centrifugal radii was only due to time evolution. Since 
a disk size of a few 100 AU is comparable to the radii 
of ac cretion disks found around classical T Tauri stars 
(e.g.. iDutrev et al"]|1996f ). our median value for R c (and 
thus outer disk radii) of 60 AU suggests that the ac- 
cretion disks found around protostars will expand, as is 
expected due to angular momentum transfer in the ac- 
cretion process. 

We note that the disk component dominates the flux in 
the 3-8 fxm range for most of the models; for about 2/3 of 
these objects the fraction of the luminosity contributed 
by the star, ri star , is either 0.1 or 0.2, i.e., 80-90% of 
the luminosity is generated by accretion. About half of 
the objects that have no dominant disk contribution in 
the 3-8 /xm range have luminosities dominated by the 
star (r} s tar=0-S>), and the other half have smaller values 
for rj s tar, but generally large inclination angles. Thus, 
the Class I objects in our sample whose luminosity is 
dominated by accretion (and which likely have large mass 
accretion rates through their disks) generate SEDs that 
are dominated by the disk in the 3-8 ^m range, except if 
the disk emission is diminished due to a high inclination 
angle. 

The reference density p\ varies from 4.0 x 10 -15 g cm -3 
to 7.0 x 10~ 14 g cm -3 , with a median value of 3.0 x 10~ 14 
g cm~ 3 , very similar to the value found by KCH93. p\ 
can be used to estimate the mass infall rate from the en- 
velope onto the disk by applying Equation[IJ for a 0.5 M Q 
object, the mass infall rates for the objects in our Class 
I sample lie between 5.3 x 10~ 7 and 9.3 x 10~ 6 M yr _1 , 
about two orders of magnitude larger than the accre- 
tion r ates of classical T Tauri stars (e.g..lGullbrin g~et alJ 
1998). This would confirm previous results that Class I 
objects are at an earlier evolutionary stage than T Tauri 
stars and stil l accreting a substantial fraction of the final 
stellar mass (jMuzerolle. Hartmann. Sz Calvetlll998l ). 

Even though the cavity semi-opening angle and, in 
some cases, the inclination angle are usually less well- 
constrained than the reference density and centrifugal ra- 
dius, we observe that our models indicate a large range of 
inclination angles, covering values from 15° to 89° (with 
a median value of 55°), and a smaller range of cavity 
semi-opening angles (0.1° to 27°, with a median of 6°). 

In Table 2] we list the luminosities of the Class I ob- 
jects in our sample and their inclination angles. Both 
the values of bolometric luminosity from the literature 
and those we derived by integrating under the SED are 
actually apparent luminosities, since envelopes are not 
spherically symmetric. Thus, we would expect the appar- 
ent bolometric luminosity to be larger than the system 
luminosity (which enters as a model parameter) at low 
inclination angles, while at high inclination angles the 
opposite should be the case. The data in Table H] shows 
that this trend is only weak for our sample of Class I ob- 
jects; in general, our model luminosities are higher than 
the apparent bolometric luminosities. This could sug- 
gest that we overestimated the luminosity of some of our 
lower-inclination sources. 

We note that most of our models do not require an 



initially flattened density distribution (rj > 1); all but 5 of 
the 22 objects we modeled can be fitted by TSC models, 
which assume an initially spherically symmetric density 
distribution. However, we cannot exclude that sheet- 
collapse model fits could also be found for the remaining 
17 Class I objects for which only TSC models have been 
considered so far. Thus, as a starting point for models, it 
seems that star-forming cores in Taurus can be described 
by the approximation of collapsing, singular isothermal 
spheres. 

In a few cases we found that an additional compo- 
nent of cold dust seems to be present, since certain mod- 
els are able to reproduce the sub-mm and mm fluxes 
measured in smaller apertures (10"-20"), but not the 
measurements at similar wavelengths derived from large 
apertures (4 0"-60"). In pa r ticula r the 1.3 mm data 
points from iMotte fc Andrei (|2001[ ). which correspond 
to fluxes integrated over a beam 60" in diameter, are 
oft en higher than the model prediction. As suggested 
bv lJavawardhana et al.l (|200fl) . cold dust of roughly con- 
stant density outside of the infall region could contribute 
to the long-wavelength emission. In these cases enve- 
lope models could aid in the distinction between emis- 
sion from infalling envelope material and emission from 
the surrounding dense cloud. 

When comparing our modeling results with previous 
efforts, we note that our results often differ considerably. 
However, in this work we have the additional constraints 
through our IRS spectra, which, when combined with 
near-IR and far-IR to mm data, allows us to determine 
the 9 parameters entering in the envelope models. In par- 
ticular the 10 fiuY silicate absorption feature allows us to 
narrow down the parameter space. A next step in mod- 
eling would involve reproducing the images of the Class 
I sources at the same time as the SEDs to better take 
into account the envelope geometry and asymmetries. 

6. CONCLUSIONS 

After analyzing the IRS spectra and SEDs of 28 Class I 
objects in the Taurus star-forming region, and generating 
envelope models for 22 of these objects, we conclude the 
following: 

• Almost all protostars we observed display ice absorp- 
tion features in their mid-IR spectra, which originate in 
the envelopes around these young stars. In particular, 
the CO2 ice feature at 15.2 /im is ubiquitous and strong 
among the Class I objects in our sample, signifying an 
origin in the cold outer envelope regions. The envelope 
is also generally responsible for a deep silicate absorp- 
tion feature at 10 /xm. The three objects with silicate 
emission features in their spectra, IRAS 04248+2612, 
04264+2433, and CoKu Tau/1 are likely Class I objects 
seen through low-density envelopes, which are in an ad- 
vanced stage of dispersal and thus considered to be more 
evolved. 

• The SEDs of Class I objects peak in the mid- to 
far-IR, where the emission is dominated by the envelope. 
Objects with flat SEDs in the infrared spectral range 
could be in transition between the Class I and II stage. 
The inclination angle of protostellar systems strongly af- 
fects the appearance of their mid-IR SED; Class I ob- 
jects seen pole-on could appear as Class II objects with 
additional long-wavelength excess emission, while edge- 
on Class I objects have SED shapes similar to Class 
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objects, but with less flux in the sub-mm and mm wave- 
length region. 

• The IRS spectra are instrumental in constraining en- 
velope model parameters, in particular the reference den- 
sity pi (which is linked to the mass infall rate) and cen- 
trifugal radius R c . The depth of the silicate feature and 
the shape of the spectrum between 12 and 30 /im provide 
constraints which can be further tightened by model fits 
to near-IR, far-IR, sub-mm, and mm observations. The 
median values for p\ and R c for the Class I objects we 
modeled are 3 x 10 -14 g cm -3 and 60 AU, respectively, 
with a relatively large range for both parameters. The 
disk contribution dominates the emission in the 3-8 pm. 
range for most of the modeled Class I objects, confirm- 
ing the importance of disk accretion luminosity at this 
early evolutionary stage. Most SEDs can be reproduced 
by TSC models, thus not requiring an initially flattened 
density distribution of the collapsing cloud core; only 5 of 
the 22 objects we modeled needed sheet-collapse models 
for an adequate fit of their SEDs. 

The IRS spectra of Class I objects in Taurus display 
a variety in appearances and add valuable constraints 
to envelope models. Fits to the observed SEDs yield 
the parameters describing their physical structure, and 
thus aid in confirming and refining our understanding of 
protostellar cloud collapse and evolution. 
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TABLE 1 

Properties of Observed Class I Objects 



Name Alt. name Multiplicity 11 Spectral L(, ; References 

Type (L Q ) 

(1) (2) (3) (4) (5) (6) 



04016+2610 


L1489 IRS 


s 


K4 


3.7 


1, 2 


04108+2803 B b 


L1495 IRS 


21.2" 




0.62 


2 


04154+2823 




s 


M0.5-M4.5 


0.33 


3, 4 


04158+2805 




s 


M5-M6 


0.20 


1, 5, 6 


04166+2706 




s 




0.5 


7 


04169+2702 




s 




0.8 


2 


04181+2654 A d 




31.3" 




0.26 


2 


D/L1 SI 4-9fi^zl R d 




31.3" 




0.25 


2 






0.30" 




1.27 


8 2 


nzL9zL£4-9R1 9 f A R ^ , ^ c 


TTtT^I TR^9 
nnoi niOi 






u.ou 


q in 1 ^9 

y, i-U, r, o, z 




njiias u 




Ml 


n ^7 

U.O 1 


1 2 


04278+2253 (A,B) 




6.8" 


G8+K7 


7.2 


1 6 


04295+2251 


L1536 IRS 


s 




0.44 


2 


04302+2247 




s 




0.34 


2 


04325+2402 (A,B,C) e 


L1535 IRS 


0.24", 8.15" 




0.9 


11, 10, 12 


04361+2547 (A,B) 


TMR 1 


0.31" 




3.8 


14, 6 


04365+2535 


TMC 1A 


s 




2.1 


4 


04368+2557 (A,B) 


L1527 IRS 


0.17" 




1.6 


15, 6 


04381+2540 (A,B) 


TMC 1 


0.6" 




0.73 


13, 4 


04489+3042 




s 


M3-M4 


0.30 


5, 1, 2 


CoKu Tau/1 (A,B) 




0.24" 


K7+M 


>0.29 


9, 1, 2 


DG Tau B 




s 




>0.02 


2 


GV Tau (A,B) 


Haro 6-10 


1.3" 


K7 


6.98 


10, 1, 2 


HH 30 




s 


M0 


>0.1 


1, 18 


HL Tau 


04287+1807 


s 


K5 


6.60 


1, 2 


IC 2087 IR 


04369+2539 


s 


K4 


3.80 


1, 2 


L1551 IRS5 (A,B) 


04287+1801 


0.3" 


G-K 


28 


16, 17, 12 


LkHa 358 


CoKu Tau/2 


s 


M5.5 


0.59 


3 


References. — (lllWhite & Hillenbrand 


J2004D; (21 IKenvon & Hartmannl (119951); OIILuhman 



I 2000fi ; (4) IMvers et al.l (119871); ( 51 ILuhmaiJ (1200611: ( 61 IKenvon et al.l (1199011 : (7) lYoung et ~ 
booai: (81 IReipurth et al.l teOOOft; (91 IPadgett et all fl!999|1 : ( 101 IDuchene et al. | (|2004f1: 



iill 



Terebev et al.l 



.1 U998I) : f 1ST 



lHartmann et all dl999h: Cl 2 ) lMotte fc Andrei ll200l|Wl3 l JApai et all J2005h: (141 
1119980; (151 ILoinard et al.| J2002I) ; (16) IRodrfguez et all (119981) ; f 171 i IKenvon et ' 
IReipurth et al.l $199$) 

Note. — Column (1) gives the name of the object, column (2) other object names commonly 
found in the literature, column (3) the multiplicity of the object (see note (a) below), column (4) 
the spectral type, column (5) the luminosity of the source taken from the literature, and column 
(6) gives the references for the data listed in the previous columns. 

a « g jj means s i n gl c star; for multiple systems, the separation between the components in arcscconds 
is listed. 

b 04108+2803 B is separated by 2l"from 041 08+2803 A. 

c 4248+2612 A and B ar e separated by 0.16" llPadgett et all 1999ft : a third component was detected 
by IDuchene et al at a distance of 4.55 and a position angle of 15°. 

d 31.3" is the separation between 04181+2654 A and B. 

e 04325+2402 is a possible triple system; the sub-arsecond bina ry may just be a single source 
surrounded by a complex reflection nebula (Hartmann et al. 1999). 
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TABLE 2 

Long-Wavelength Data of Class I Objects 



Name 


A 


r lux 


Aperture size 


Reference 




w 


( J y) 


Vv nivi J 




04016+2610 


160 


46.0 


50" 


i 




350 


12.48 


9" 


2 




450 


4.23 


40" 


4 




800 


0.58 


16'.'8 


5 




850 


0.59 


40" 


4 




1100 


0.18 


18'.'5 


5 




1300 


0.15 


60" 


6 


04108+2803 B 


450 


1.13 


40" 


4 




800 


0.085 


16'.'8 


5 




850 


0.17 


40" 


4 




1100 


< 0.1 


18"5 


5 




1300 


0.04 


60" 


6 


04154+2823 


350 


0.44 


9" 


2 




450 


0.495 


9" 


2 




800 


<0.1 


16'.'8 


5 




850 


0.14 


15" 


2 




1100 


<0.1 


18'.'5 


5 


04158+2805 


880 


0.067 


~ 5" 


12 




1300 


0.11 


60" 


6 


04166+2706 


350 


6.94 


9" 


2 




450 


4.2 


40" 


4 




800 


0.65 


16'.'8 


3 




850 


1.08 


40" 


4 




1100 


0.30 


18"5 


3 




1300 


0.80 


60" 


6 


04169+2702 


350 


7.34 


9" 


2 




450 


6.09 


40" 


4 




800 


0.75 


16" 8 


5 




850 


1.14 


40" 


4 




1100 


0.28 


18'.'5 


5 




1300 


0.73 


60" 


6 


04181+2654 A 


1300 


0.23 


60" 


6 


04181+2654 B 










04239+2436 (A,B) 


350 


1.14 


9" 


2 




450 


<0.66 


9" 


2 




800 


0.33 


16'/8 


5 




850 


0.21 


15" 


2 




1100 


0.11 


18'/5 


5 




1300 


0.17 


60" 


6 


04248+2612 (A,B,C) 


350 


1.18 


9" 


2 




450 


2.96 


40" 


4 




800 


0.25 


16'/8 


5 




850 


0.56 


40" 


4 




1100 


0.10 


18'.'5 


5 




1300 


0.45 


60" 


6 


04264+2433 


450 


0.63 


40" 


4 




850 


0.13 


40" 


4 




1300 


0.03 


60" 


6 


04278+2253 (A,B) 


450 


<0.69 


9" 


2 


850 


0.04 


15" 


2 


04295+2251 


350 


1.34 


9" 


2 




450 


2.66 


40" 


4 




800 


0.24 


16'.'8 


5 




850 


0.42 


40" 


4 




1100 


0.09 


18'.'5 


5 




1300 


0.115 


60" 


6 


04302+2247 


350 


2.87 


9" 


2 




450 


2.09 


40" 


4 




800 


0.34 


16'/8 


5 




850 


0.57 


40" 


4 




1100 


0.15 


18'/5 


5 




1300 


0.18 


60" 


6 


04325+2402 (A,B,C) 


160 


38.0 


50" 


1 




450 


0.61 


9" 


2 




800 


0.30 


16'/8 


5 




850 


0.19 


15" 


2 




1100 


0.07 


18"5 


5 




1300 


0.52 


60" 


6 


04361+2547 (A,B) 


450 


2.35 


40" 


4 


800 


0.63 


16'.'8 


5 




850 


0.64 


40" 


4 




1100 


0.19 


18'/5 


5 




1300 


0.44 


60" 


6 
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TABLE 2 — Continued 



Name 


A 


Flux 


Aperture size Reference 






(Jy) 


(FWHM) 



04365+2535 


350 


20.6 


45" 


10 




450 


12.7 


45" 


10 




800 


1.01 


16'.'8 


5 




850 


1.80 


45" 


10 




1100 


0.44 


18'.'5 


5 




1300 


0.45 


60" 


6 


04368+2557 (A,B) 


160 


69.0 


50" 


1 


350 


12.0 


45" 


1 




450 


2.85 


9" 


2 




800 


1.52 


16"8 


5 




850 


0.895 


15" 


2 




1100 


0.48 


18"5 


5 




1300 


1.50 


60" 


6 


04381+2540 (A,B) 


450 


2.82 


40" 


i 




800 


0.29 


16V8 


5 




850 


0.56 


40" 


4 




1100 


0.12 


18"5 


5 




1300 


0.30 


60" 


6 


04489+3042 


1300 


> 0.015 


60" 


6 


CoKu Tau/1 (A,B) 


450 


<0.52 


9" 


2 




850 


0.035 


15" 


2 




1300 


< 0.012 


11" 


9 


DG Tau B 


1300 


0.31 




11 


GV Tau (A,B) 


350 


1.68 


9" 


2 


450 


1.81 


9" 


2 




800 


0.57 


16" 


8 




850 


0.28 


15" 


2 




1100 


0.18 


19" 


8 




1300 


0.20 


60" 


6 


HH 30 


1300 


0.035 


60" 


6 


HL Tau 


350 


26.4 


45" 


10 




450 


16.8 


45" 


10 




800 


2.58 


15"8 


7 




850 


2.97 


45" 


10 




1100 


1.11 


18^4 


7 




1300 


1.20 


60" 


6 


IC 2087 IR 


450 


1.365 


9" 


2 




850 


0.50 


15" 


2 


L1551 IRS5 (A,B) 


350 


164 


45" 


10 


450 


94 


45" 


10 




800 


8.05 


16^8 


5 




850 


12.1 


45" 


10 




1100 


2.77 


18"5 


5 




1300 


3.40 


60" 


6 


LkHa 358 


1300 


0.032 


11" 


9 



References. (1) ILadd et alj U991f: (2) [ Andrews h Williams] 120051); ( 3) IBarsonv fc Kenvonl IU992T ): (4) lYoung et al.l ||2003J1: (5 ) 

IMoriartv-Schieven et al I ||1994 ): (6 HMotte &: Andrei ihuOlf ): (71 [Adams et alj dl990ft: (81 | Chandler et al.l (Il998f l: (9) IQsterloh fc Beckwithl lll995f l: 
(10) IChandler fc Riche? II2000H : fill IPadgett et al.l IU999T) : ( 12) Undrews fc Williamd 112007^ 
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TABLE 3 
Model Fits for Class I Objects 



Name 


L 


Pi 


R c 




^max 


V 


Vstar 


C0 2 ice e i 




(Lq) 


(g cm- 3 ) 


(AU) 


(F^star ) 


(AU) 






abundance 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) (10) (11) 



04016+2610 


4.5 


4.5 x 10" 


14 


100 


1 


6000 


1.0 


0.1 


1.5 x 


nr 


-4 


5° 


40° 


04108+2803 B 


0.7 


1.5 x 10~ 


14 


40 


3 


6000 


TSC 


0.3 


1.0 x 


10" 


-4 


10° 


40° 


04154+2823 


0.35 


7.0 x 10~ 


15 


10 


1 


5000 


TSC 


0.1 


1.0 x 


io- 


-4 


5° 


20° 


04158+2805 


0.3 


2.0 x 10" 


14 


60 


3 


5000 


TSC 


0.7 


4.0 x 


io- 


-5 


5° 


30° 


04166+2706 


0.6 


4.5 x 10- 


14 


300 


5 


10000 


TSC 


0.2 


2.0 x 


io- 


-5 


6° 


85° 


04169+2702 


1.5 


3.2 x 10- 


14 


100 


1 


10000 


1.0 


0.8 


6.0 x 


io- 


-5 


0.1° 


75° 


04181+2654 A 


0.7 


2.0 x 10- 


14 


50 


5 


10000 


TSC 


0.2 


7.0 x 


10" 


-5 


7° 


20° 


04239+2436 


1.5 


1.3 xlO- 
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Note. — Column (1) gives the name of the object, column (2) the total luminosity (L s t ar + ^dis/c) of the system, column (3) the reference 
density pi, column (4) the centrifugal radius, column (5) the inner disk radius (1 ~R s tar — 2 Rq), column (6) the outer radius of the envelope, 
column (7) gives the flattening parameter r/, column (8) the fraction of the luminosity arising from the star, column (9) the CO2 ice abundance, 
column (10) the semi-opening angle of the cavity, and column (11) the inlincation angle. 

Typical uncertainties for our model parameters are as follows: L ± 10%, p\ ± 15%, R c ± 20%, R-diskmin =t 10%, R mQE ± 25%, rj ± 0.25 (if not 
TSC model), r) staT ± 0.05, CO2 ice abundance ± 10%, 9 ± 2°, i ± 5°. Note that most uncertainties arc correlated; e.g., a decrease in pi should 
be accompanied by a decrease in R c and/or increase in i to still result in a comparable fit. 
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TABLE 4 

Comparison of Luminosities and Inclination Angles 



Name 


OOl \ \y / 
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0.35 






20 


04158+2805 


0.20 


3 


0.2 


0.3 






30 


04166+2706 


0.5 (0.40) 


4(3) 


0.4 


0.6 






85 


04169+2702 


0.8 (1.4) 


1 (3) 


1.2 


1.5 


60 


10 


75 


04181+2654 A 


0.26 (0.70) 


1 (3) 


0.7 


0.7 






20 


04181+2654 B 


0.25 


1 


0.3 










04239+2436 (A,B) 


1.27 


1 


1.5 


1.5 






15 


04248+2612 (A,B,C) 


0.36 


1 


0.2 


0.4 


78 


8 


70 


04264+2433 


0.37 


1 


0.5 


0.7 






87 


04278+2253 (A,B) 


7.2 


3 


5.7 










04295+2251 


0.44 (0.64) 


1 (2) 


0.5 


0.8 






70 


04302+2247 


0.34 


1 


0.3 


1.0 


90 


8 


89 


04325+2402 (A,B,C) 


0.9 (0.70) 


5(1) 


1.2 


0.9 


60 
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3.8 (2.90) 


3(1) 
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60 


7 
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2.5 
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7 


30 


04368+2557 (A,B) 
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1.8 


75 


7 
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0.73 


2 


0.7 


1.0 


55 


7 
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0.30 


1 


0.3 


0.3 






20 


CoKu Tau/1 (A,B) 


>0.29 


1 


1.0 


1.1 


~90 


9 


80 


DG Tau B 


>0.02 


1 


1.8 


2.5 


~90 


9 
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GV Tau (A,B) 


6.98 


1 


9.0 




30 


7 




HH 30 


>0.1 


6 


0.02 




~90 


12 




HL Tau 


6.60 (7.1) 


1 (6) 


6.9 


8.0 


67 


11 


15 


IC 2087 IR 


3.80 




4.9 










L1551 IRS5 (A,B) 


28 (21.90) 


5(1) 


23.2 


25.0 


65 (56) 


7(8) 


45 


LkHa 358 


0.59 


13 


0.3 











Ref eren ces. — (DlKenyon & Hartmannl J1995I); (21 iMyers et alii 198 71): QllKenyon et all dl990l) (4)|Young et all d2003|); (5)|Motte & Andrei 
J2u0l[); (6)lReipurth et al.||l993D; (7) iHogerheiide et al.l J1998D; f8nLucas Roche! dl997l) ; (9) iPadgett et al.l J1999I7 : flOl lOhashi et all dl997b|) ; 
qil lClose et all d!997l) ; Cl2T lBurrows et al.l J1996D ; fl3) |Luhman| d2000|) 

Note. — Column (1) gives the name of the object, columns (2) and (3) the bolometric luminosity of the source taken from the literature and 
its reference, respectively, column (4) gives the bolometric luminosity measured by integrating under the SED, column (5) the luminosity used 
in the models, columns (6) and (7) give the inclination angle taken from the literature and its reference, respectively, and column (8) gives the 
inclination angle used in the models. 
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Fig. 1. — The IRS SL, LL, SH, and LH slit po s itions (only one of the two nod positions is shown) superposed on the NICMOS J-band 
(F110W) image of DG Tau B from Padgett et al. ( 199EJ); note that for this object the LL module was not used and is shown for illustrative 
purposes only. The slits, with their different widths and orientations, cover different parts of the object. 
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Fig. 2. — IRS Spectra of Class I Objects in Taurus. 
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Fig. 2. — continued. 
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Fig. 2. — continued. 
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Fig. 2. — continued. 
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Fig. 2. — continued. 
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Fig. 3. — SED plots of the Class I objects in our sample, ordered alphabetically by their name. Optical to mid-IR, ground-based 
photometry is shown as filled diamonds, the 2MASS J, H, and K s fluxes as filled circles, the IRAC 3.6, 4.5, 5.8, and 8.0 fim fluxes as open 
circles, the IRAS 12, 25, and 60 fim fluxes as open diamonds or open, upside down triangles, if upper limit, and sub-mm and mm fluxes, 
where available, as filled squares. The IRS spectrum is also shown. No corrections for reddening have been applied. 
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Fig. 3. — continued. 
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Fig. 3. — continued. 
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Fig. 4. — A sketch of an envelope model as described in S I4.1I 
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Fig. 5. — IRS spectrum and photometric data of 04016+2610, and an envelope model fit with L=4.5 Lq, pi = 4.5 X 10 g cm - , 
R c =100 AU, n=l.Q, rj s tar=0.l, 8=5°, i=40°, an inner disk radius of 1 stellar radius, and an outer envelope radius of 6000 AU. The CO2 
ice abundance set to 1.5 X 10~ 4 . The gray lines represent the different model components: the envelope (long-dashed line), the star (dotted 
line) and the disk (dash-dotted line), both extinguished by the envelope, and the sum of all components (solid gray line). 
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Fig. 6.— IRS spectra of both components of the 21" binary IRAS 04108+2803, plotted on the same scale. While observing 04108+2803A, 
the bright B component partly entered the LL slit, and therefore the spectrum of component A beyond 14 (im contains an increasing amount 
of emission from component B; the spectrum of A is dominated by emission from B beyond about 20 /im. 
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Fig. 7.— 
R c =40 AU, 
of 6000 AU. 
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IRS spectrum and photometric data of 04108+2803B, and an envelope model fit with L=0.7 Lq, pi = 1.5 X 10~ 
initial TSC density distribution, Jjstar=0.3, 9=10°, i=40°, an inner disk radius of 3 stellar radii, and an outer envelope radius 
The CO2 ice abundance was set to 1.0 X 10~ 4 . The gray lines represent the model components as in Figure[5] 
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Fig. 8. — IRS spectrum and photometric data of 04154+2823, and an envelope model fit with L=0.35 L©, pi = 7.0 X 10 15 g cm 3 , 
_R C =10 AU, initial TSC density distribution, fjstar — 0.1, 9 — 5°, i — 20°, an inner disk radius of 1 stellar radius, and an outer envelope radius 
of 5000 AU. The CO2 ice abundance was set to 1.0 x 10 -4 . The gray lines represent the model components as in Figure[5] 
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Fig. 9.— 
R c =60 AU, 
of 5000 AU. 
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L=0.3 L sun , p,=2.0e-14 g/cm 3 , R c =60 AU, TSC, T) stQr =0.7, 6 = 5°, i = 30° 




1000 



IRS spectrum and photometric data of 04158+2805, and an envelope model fit with L=0.3 Lq, pi = 2.0 X 10 — 14 g cm -3 , 
initial TSC density distribution, r) a tar=0.7, 8=5°, i=30°, an inner disk radius of 3 stellar radii, and an outer envelope radius 
The CO2 ice abundance was set to 4.0 X 10 — 5 . The gray lines represent the model components as in Figure [5] 
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Fig. 10. — IRS spectrum and photometric data of 04166+2706, and an envelope model fit with L=0.6 Lq, pi = 4.5 X 10 — 14 g cm -3 , 
-R c =300 AU, initial TSC density distribution, ?j 3 t ar =0.2, 0=6°, i=85°, an inner disk radius of 5 stellar radii, and an outer envelope radius 
of 10000 AU. The CO2 ice abundance was set to 2.0 x 10 -5 . The gray lines represent the model components as in Figure [5] 
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L=1.5 L sun , p,=3. 2e-14 g/cm 3 , R c =100 AU, 77 = 1.0, 7) stor =0.8, 6 = 0.1°, i = 75° - 
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Fig. 11. — IRS spectrum and photometric data of 04169+2702, and an envelope model fit with L=1.5 L0, pi = 3.2 X 10 — 14 g cm -3 , 
-R c =100 AU, 77=1.0, 77 s t ar =0.8, 9=0.1°, i=75°, an inner disk radius of 1 stellar radius, and an outer envelope radius of 10000 AU. The CO2 
ice abundance was set to 6.0 x 10 -5 . The gray lines represent the model components as in Figure [5] 
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Fig. 12. — IRS spectrum and photometric data of 04181+2654A, and an envelope model fit with L=0.7 Lq, pi = 2.0 X 10 -14 g cm -3 , 
R c =50 AU, initial TSC density distribution, r/ s t ar =0.2, 8=7° , i=20°, an inner disk radius of 5 stellar radii, and an outer envelope radius 
of 10000 AU. The CO2 ice abundance was set to 7.0 x 10 -5 . The gray lines represent the model components as in Figure [5] 
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Fig. 13. — IRS spectrum and photometric data of 04239+2436, and an envelope model fit with L=1.5 L©, pi = 1.3 X 10 14 g cm 3 , 
_R C =10 AU, initial TSC density distribution, ^stai — 0.1, 9 — 5°, i — 15°, an inner disk radius of 1 stellar radius, and an outer envelope radius 
of 5000 AU. The CO2 ice abundance was set to 1.0 X 10~ 4 . The gray lines represent the model components as in Figure[5] 
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Fig. 14. — IRS spectrum and photometric data of 04248+2612, and an envelope model fit with L=0.4 Lq, pi = 4.0 X 10 5 g cm -3 , 
ij c =30 AU, initial TSC density distribution, rj B tar=0.3, 8=15° , i=70°, an inner disk radius of 5 stellar radii, and an outer envelope radius 
of 10000 AU. The CO2 ice abundance was set to 5.0 x 10 -5 . The gray lines represent the model components as in Figure [5] 
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L=0.7 L sun , p,=5.0e-15 g/cm 3 , R c =30 AU, TSC, 77 stor =0.8, 6=13°, i = 87° 




Fig. 15. — IRS spectrum and photometric data of 04264+2433, and an envelope model fit with L=0.7 Lq, pi = 5.0 X 10 — 15 g cm -3 , 
ij c =30 AU, initial TSC density distribution, r) s tar=0.8, 6=13°, i=87°, an inner disk radius of 5 stellar radii, and an outer envelope radius 
of 10000 AU. The CO2 ice abundance was set to 5.0 x 10 -5 . The gray lines represent the model components as in Figure [5] 
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Fig. 16. — IRS spectrum and photometric data of 04295+2251, and an envelope model fit with L=0.8 Lq, pi = 8.0 X 10 15 g cm 3 , 
_R C =20 AU, initial TSC density distribution, flstai — 0.2, 6 — 5° j i — 70°, an inner disk radius of 5 stellar radii, and an outer envelope radius 
of 10000 AU. The CO2 ice abundance was set to 1.0 X 10 -4 . The gray lines represent the model components as in Figure [5] 
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Fig. 17. — IRS spectrum and photometric data of 04302+2247, and an envelope model fit with L=1.0 Lq, pi = 3.0 X 10 -14 g cm -3 , 
-R c =300 AU, initial TSC density distribution, r) s t ar =0.3, 9=22°, i=89°, an inner disk radius of 1 stellar radius, and an outer envelope 
radius of 10000 AU. The CO2 ice abundance was set to 4.0 X 10 — . The gray lines represent the model components as in Figure |51 




Fig. 18.— The IRS SL and LL slit positions superposed on the NICMOS K-band (F205W) image of IRAS 04325+2402 by Hartmann ct al. 
(1999). The SL slit (3"6 wide) runs from the lower left to upper right, while the wider LL slit extends from the upper left to the lower 
right. The full slit lengths are larger than the size of the image. 
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Fig. 19. — IRS spectrum and photometric data of 04325+2402, and an envelope model fit with L=0.9 Lq, pi = 3.0 X 10 -14 g cm -3 , 
-R c =100 AU, 77=1.0, ristar=0.8, 0=15°, i=80°, an inner disk radius of 3 stellar radii, and an outer envelope radius of 5000 AU. The CO2 
ice abundance was set to 4.0 x 10 -5 . The gray lines represent the model components as in Figure [5] 
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Fig. 20. — IRS spectrum and photometric data of 04361+2547, and an envelope model fit with L=4.0 Lq, pi = 2.0 X 10 -14 g cm -3 , 
-R c =100 AU, T] =1.5, r) B tar=0.9, 9=15°, i=80°, an inner disk radius of 2 stellar radii, and an outer envelope radius of 10000 AU. The CO2 
ice abundance was set to 7.0 x 10 -5 . The gray lines represent the model components as in Figure [5] 
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Fig. 21. — IRS spectrum and photometric data of 04365+2535, and an envelope model fit with L=2.5 Lq, pi = 4.5 X 10 — 14 g cm -3 , 
-R c =50 AU, initial TSC density distribution, rj B tar=0.1, 6=5°, i=30°, an inner disk radius of 1 stellar radius, and an outer envelope radius 
of 10000 AU. The CO2 ice abundance was set to 7.0 x 10 -5 . The gray lines represent the model components as in Figure [5] 
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Fig. 22. — IRS spectrum and photometric data of 04368+2557, and an envelope model fit with L=1.8 Lq, pi = 4.0 X 10 — " g cm" 
-R c =200 AU, initial TSC density distribution, r? s tar=0.5, 9=27°, i=89°, an inner disk radius of 1 stellar radius, and an outer envelope 
radius of 10000 AU. The CO2 ice abundance was set to 1.0 X 10 — 4 . The gray lines represent the model components as in Figure |51 
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Fig. 23. — IRS spectrum and photometric data of 04381+2540, and an envelope model fit with L=1.0 Lq, pi = 3.0 X 10 — 14 g cm -3 , 
R c =70 AU, initial TSC density distribution, r) s tar=0.1, 6=10°, i=40°, an inner disk radius of 5 stellar radii, and an outer envelope radius 
of 10000 AU. The CO2 ice abundance was set to 1.2 X 10 — 4 . The gray lines represent the model components as in Figure[5] 
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Fig. 24. — IRS spectrum and photometric data of 04489+3042, and an envelope model fit with L=0.3 Lq, pi = 1.0 X 10~ 14 g cm -3 , 
R c =15 AU, initial TSC density distribution, r) B tar=0.7, 8=1°, i=20°, an inner disk radius of 1 stellar radius, and an outer envelope radius 
of 1000 AU. The CO2 ice abundance was set to 4.0 X 10 . The gray lines represent the model components as in Figure [5] 
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Fig. 25. — IRS spectrum and photometric data of CoKu Tau/1, and an envelope model fit with L=l.l L©, pi = 5.0 X 10 15 g cm 3 , 
R c =40 AU, initial TSC density distribution, Vstai — 0.8, 6 — 5°, i — 80°, an inner disk radius of 7 stellar radii, and an outer envelope radius 
of 5000 AU. The CO2 ice abundance was set to 5.0 X 10 -5 . The gray lines represent the model components as in Figure [5] 
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Fig. 26. — IRS spectrum and photometric data of DG Tau B, and an envelope model fit with L=2.5 Lq, pi = 3.5 X 10 — 14 g cm -3 , 
ij c =60 AU, initial TSC density distribution, ?7star=0.1, 0=10°, i=55°, an inner disk radius of 1 stellar radius, and an outer envelope radius 
of 10000 AU. The CO2 ice abundance was set to 3.0 x 10 -5 . The gray lines represent the model components as in Figure [5] 
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Fig. 27. — IRS spectrum and photometric data of HL Tau, and an envelope model fit with L=8.0 Lq, p\ = 4.5 X 10~ 14 g cm -3 , 
-R c =100 AU, r;=1.0, r] s tar=0.2, 8=1°, i=15°, an inner disk radius of 2 stellar radii, and an outer envelope radius of 10000 AU. The CO2 
ice abundance was set to 4.0 x 10 -5 . The gray lines represent the model components as in Figure [5] 
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Fig. 28. — IRS spectrum and photometric data of L1551 IRS 5, and an envelope model fit with L=25.0 Lq, pi = 7.0 X 10 — 14 g cm -3 , 
-R c =100 AU, initial TSC density distribution, r] s tar=0.3, 8=5°, i=45°, an inner disk radius of 5 stellar radii, and an outer envelope radius 
of 10000 AU. The CO2 ice abundance was set to 3.0 x 10 -4 . The gray lines represent the model components as in Figure [5] 



4G 



Furlan et al. 



CD 
3 



O 
CD 



o 
o 



csl 

Lf) 



2.0 



.5 



+ 



A 



0.5 * 

m 



i 



o.o gl 



A 
B 

A A 



♦ Taurus Class I objects 
Taurus Class l/ll objects 
A embedded, massive YSOs 

(Gibb et al. 2004) 
□ low- and intermediate-mass YSOs 

(Alexander et al. 2003) 



peak 



(10 fxru silicate feature) 



Fig. 29. — Peak optical depth of the CO2 ice feature versus the peak optical depth of the 10 /im silicate absorption feature. The diamonds 
represent the Cla ss I objects presented in thi s paper, while the triangles and squares represent generally more massive YSOs studied by 
Gibb ct al. (2004) and Alexander ct al. (2003), respectively. The data points at a peak silicate optical depth of are objects with either a 
silicate emission feature or a feature that shows emission and absorption characteristics. The thick, gray dashed line is drawn to guide the 
eye to separate the regimes of envelopes from that of ambient molecular cloud material (see text for details). 



